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ABSTRACT 

A  computer  controlled  data  acquisition  system  utilizing 
hotwire  anemometry  has  been  designed,  built,  and  installed 
in  the  Naval  Postgraduate  School  Low  Speed  Wind  Tunnel 
Facility.  All  relevant  wind  tunnel  data  is  obtained  by  the 
use  of  a  computer  guided  Hewlett  Packard  data  acquisition 
system  and  a  Para-Scientific  Pressure  Computer.  Two  computer 
programs  were  written  to  coordinate  hotwire  system 
calibration  with  hotwire  experimentation.  An  experiment, 
Wake  Velocity  Profile  Analyisis  and  Drag  Coefficient 
Measurement  of  an  Airfoil,  was  used  as  a  vehicle  to  test  the 
system.  The  final  output  of  the  data  acquisition  system, 
including  graphical  information,  compared  favorably  with 
previous  results  from  an  older  data  acquisition  system 
already  in  use.  Drag  coefficient  output  compared  very 
favorably  to  data  supplied  by  the  National  Advisory 
Committee  on  Aeronautics. 
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I.  INTRODUCTION 

It  was  the  goal  of  this  thesis  project  to  design,  build 
and  implement  a  computer  operated  data  acquisition  system 
for  the  NPS  wind  tunnel  utilizing  hotwire  anemometry. 
Preliminary  project  goals  were  to  interface  the  new  hotwire 
data  acquisition  system  with  the  existing  pitot-static 
system  for  on-line  computer  analysis  and  comparison. 
Additionally,  the  new  computer  program  must  produce  hardcopy 
output  of  all  graphical  and  tabular  data.  A  final  design 
requirement  was  storage  of  all  wake  velocity  profile  data 
for  future  analysis.  Versatility  was  the  foremost  design 
requirement  for  the  data  acquisition  computer  program.  The 
computer  program  is  to  have  application  as  a  classroom 
demonstration  tool  or  laboratory  introduction  to  hotwire 
anemometry.  The  system  can  also  be  used  in  future  research 
of  aerodynamic  models  or  airfoil  sections. 

To  enhance  adaptability,  the  computer  program  was 
created  in  two  parts.  The  first  part,  program  CALIB, 
performs  system  set  up  and  calibrates  the  hotwire.  As  a  test 
case  for  the  part  two  application  program,  an  experiment  was 
conducted  where  the  wake  velocity  profile  of  a  NACA  66(215)- 
216  airfoil  was  analyzed. 
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Both  parts  of  the  application  program  guide  the  user 
through  each  step  of  the  experiment.  Raw  data,  reduced  data, 
and  graphical  outputs  are  all  tabulated  on  screen  for 
immediate  examination.  Hard  copy  output  of  all  on  screen 
information  is  also  available.  The  second  application 
program,  PART 2 ,  provides  wake  velocity  profiles  of  the 
airfoil  using  hotwire  anemometry  and  a  pitot-static  system 
for  comparison.  Drag  coefficients  are  also  calculated  for 
the  airfoil  as  a  function  of  air  flow  velocity  and  angle  of 
incidence. 


II.  THEORY  AND  INSTRUMENTATION 

A.   HOT  WIRE  ANEMOMETRY 

A  thermal  anemometer  is  used  to  measure  fluid  velocity 
by  responding  to  changes  in  heat  transfer  from  an 
electrically  heated  probe  immersed  in  a  fluid  flow.  The 
amount  of  electrical  energy  dissipated  from  the  sensor  is  a 
measure  of  the  cooling  effect  of  the  fluid  flowing  past  the 
heated  sensor.  Heat  loss  from  the  hotwire  probe  is  dependent 
on  the  wire  operating  temperature,  geometry  of  the  wire, 
thermal  properties  of  the  wire,  and  physical  properties  of 
the  fluid  flow  (ie,  air  temperature,  velocity,  pressure,  and 
density) . 

Air  flow  in  the  NPS  wind  tunnel  can  be  considered  low 
speed,  steady,  incompressible  flow.  Under  these  conditions  a 
reasonable  assumption  can  be  made  that  the  fluid 
temperature,  density  and  pressure  are  non-varying.  With  this 
assumption  heat  loss  from  the  wire  sensor  can  be  interpreted 
as  a  direct  measure  of  the  fluid  velocity. 

The  sensing  element,  also  called  the  hotwire  sensor  or 
probe,  is  a  fine  metal  wire  (4-5  microns  in  diameter) 
suspended  between  two  metal  supports.  A  cylindrical  hot  film 
sensor  consists  of  a  wire  shaped  electrical  insulator  called 
a  substrate  with  a  conducting  film  deposited  on  the  outside 
surface.  Its  outward  appearance  is  very  similar  to  that  of  a 
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hotwire  (see  Figure  1) ,  the  main  difference  being  diameter 
(about  50  microns  for  a  hot  film  sensor  as  compared  to  4-5 
microns  for  a  hotwire) . 

The  sensor  is  mounted  on  a  long  thin  probe  in  an  effort 
to  be  as  unobtrusive  to  the  ambient  flow  as  possible. 
Electrical  current  to  the  sensor  is  supplied  by  the 
anemometer,  consisting  of  a  Wheatstone  Bridge  and  a  control 
amplifier  (see  Figure  2) .  The  other  resistances  in  the 
bridge  (R2  and  R3)  are  fixed. 

The  wheatstone  bridge  and  control  amplifier  work  on  a 
voltage  balancing  principle.  Air  flow  in  the  wind  tunnel 
will  cool  the  sensor  resulting  in  a  decrease  in  its 
resistance.  The  decreased  resistance  of  the  probe  lowers  the 
voltage  across  the  wheatstone  bridge  which  in  turn  lowers 
the  input  to  the  control  amplifier.  The  control  amplifier  is 
phase  adjusted;  a  lower  input  voltage  causes  a  current 
increase  from  the  amplifier  to  the  bridge.  The  increased 
current  heats  up  the  sensor  which  also  increases  its 
resistance.  When  the  probe's  resistance  is  back  to  its  pre- 
set value,  voltages  across  the  bridge  are  balanced  and  the 
control  amplifier  no  longer  alters  current  flow  through  the 
bridge.  If  the  amplifier  gain  is  set  correctly,  voltage 
inputs  will  not  deviate  far  from  the  balanced  condition. 
Therefore  any  change  in  sensor  resistance  due  to  fluid  flow 
fluctuations  will  be  corrected  immediately. 
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TSI  Model  1210-20  Hot  Film  Sensor 
Figure  1 
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Schematic  of  a  Constant  Temperature  Hotwire  System 

.Figure  2 
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The  control  amplifier  is  able  to  maintain  a  voltage 
balance  across  the  hotwire  probe  and  wheatstone  bridge  at  a 
rate  of  almost  2  00,000  cycles  per  second.  Thus,  hotwire 
anemometry  is  ideally  suited  for  rapidly  changing  or 
turbulent  flow  regimes.  Quick  response  to  velocity 
fluctuations  also  provides  an  accurate,  repeatable  output 
when  steady  flow  conditions  prevail.  The  IFA  100  Thermal 
Anemometer  used  in  this  thesis  is  shown  in  Figure  3. 

B.   HEAT  TRANSFER  IN  A  HOTWIRE  SENSOR 

Electrical  heat  input  to  a  hotwire  probe  is  dissipated 
in  various  modes  of  heat  transfer.  Radiation,  buoyant 
convection,  conduction  to  the  sensor's  end  supports  and 
forced  convection  all  contribute  to  sensor  cooling  in  the 
flow  field.  Radiative  heat  loss  from  a  typical  hotwire 
sensor  is  less  than  0.1  percent  of  the  heat  input  to  the 
probe.  In  very  low  density  gaseous  flows  this  source  of  heat 
loss  may  be  important.  Radiative  heat  losses  were  considered 
negligible  in  this  thesis. 

Buoyant  convection  is  important  in  heat  loss 
calculations  only  at  very  low  wind  tunnel  speeds.  Buoyant 
convection  can  be  neglected  at  fluid  speeds  greater  than  5 
centimeters/second  (about  one  mile  per  hour) .  The  primary 
contributors  to  heat  transfer  from  a  hotwire  sensor  are 
conduction  to  the  end  supports  and  forced  convection  by  the 
fluid  flow. 
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The  supports  that  hold  the  hotwire  sensor  to  the  probe 
legs  are  much  thicker  than  the  wire  sensor.  This  serves  a 
dual  purpose:  (1)  The  supports  contribute  to  the  sensor's 
strength;  (2)  The  supports  should  not  be  heated  appreciably 
by  an  electrical  current  flowing  through  the  wire  sensor. 
When  the  sensor  is  operating  in  a  flow  field  the  temperature 
along  the  wire  is  greatest  in  the  center.  The  temperature  at 
the  center  is  closest  to  that  of  an  infinite  wire,  not 
affected  by  the  supports.  Nearer  to  the  supports  the  wire 
temperature  decreases  until  it  is  very  close  to  the 
temperature  of  the  ambient  flow  field.  Heat  transfer  from 
the  sensor  to  the  supports  by  conduction  is  a  loss  and  a 
potential  source  of  error.  " 

It  has  been  experimentally  determined  that 
approximately  17  percent  of  the  heat  input  to  a  typical 
hotwire  sensor  is  lost  to  the  supports  through  conduction. 
The  remaining  heat  input  dissipates  through  the  process  of 
forced  convection.  These  numbers  represent  large 
experimental  errors  if  the  heat  loss  modes  are  not  accounted 
for.  Heat  conduction  to  the  end  supports  is . the  best 
argument  for  calibrating  the  hotwire  system  with  the  wind 
tunnel  on  (under  operating  conditions)  instead  of  using 
universal  calibration  constants.  The  calibration  curve,  if 
created  under  operating  conditions,  will  properly  account 
for  conduction  losses. 
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C.   HP  DATA  ACQUISITION  SYSTEM 

The  HP  Data  Acquisition  System  consists  of  a  combination 
of  hardware  and  software  that  allows  the  IBM   PC-AT  computer 
to  act  as  a  fully  automated  instrumentation  system.  The 
individual  HP  instruments  used  in  this  thesis  include  the 
Relay  Multiplexer,  Digital  Multimeter,  Relay  Actuator,  and 
Digital  Input  Port  (see  Figure  4) .  Each  of  the  instruments 
can  be  operated  manually  at  the  computer  screen  by  means  of 
the  mouse  control.  This" is  made  possible  by  a  software 
program  called  the  PC  Instruments  Soft  Front  Panel.  In  this 
mode  the  instruments'  operating  controls  and  functions,  as 
well  as  digital  displays  are  relayed  to  the  computer  screen. 
Instead  of  user  interaction  with  the  instruments'  controls 
and  indicators,  each  unit  is  manipulated  by  the  computer 
mouse  controller. 

In  addition  to  the. manual  mode,  the  HP  data  acquisition 
system  can  operate  in  the  programmed  (basic)  mode.  This 
method  is  employed  in  this  thesis.  A  program  is  written  in 
the  Advanced  Basic  language  (BASICA)  that  calls  up  each 
instrument  and  its  function  as  necessary.  All  data 
acquisition  related  to  the  hotwire  wind  tunnel  experiment  is 
accomplished  by  the  use  of  a  BASICA  program.  User 
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IFA  100  Thermal  Anemometer 
Figure  3 


HP  Data  Acquisition  System 
Figure  4 
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interaction  in  the  experiment  is  limited  to  those  tasks  not 
controlled  by  the  HP  data  acquisition  system. 

D.   INSTRUMENTATION 

1.  Relay  Multiplexer 

The  Relay  Multiplexer  is  capable  of  switching 
between  eight  input  channels  of  voltage  data.  Each  of  the 
eight  inputs  consists  of  two  relay  contacts  (double  pole- 
single  throw) .  The  system  provides  one  common  output  channel 
which  in  turn  can  be  read  by  the  Digital  Multimeter.  Multple 
relay  closures  cannot  occur;  only  one  input  can  be  switched 
to  the  common  output  at  a  time.  The  Relay  Multiplexer 
utilizes  break-bef ore-make  input  switching.  The  presently 
selected  input  channel's  signal  is  opened  and  interrupted 
before  a  new  input  channel  is  connected  to  the  common 
output.  Figure  10  in  Chapter  III  shows  all  data  signal  paths 
routed  to  the  Relay  Multiplexer. 

2 .  Digital  Multimeter 

The  Digital  Multimeter  (DMM)  can  measure  +  or  -  DC 
voltages,  AC  voltages,  and  ohms.  In  this  thesis  only  the  DC 
voltage  measurement  function  is  used.  Direct  current 
voltages  can  be  read  on  four  ranges  varying  from  0.2  V  to 
2  00  V.  An  auto  range  feature  also  exists  where  the  DMM 
selects  the  optimal  range  for  the  signal  that  is  being 
measured. 
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The  DMM  automatically  converts  input  analog  voltage 
signals  into  a  digital  (or  binary)  form  which  can  be  read  by 
the  computer.  When  the  DMM  sends  its  measurement  to  the 
computer  its  value  is  assigned  to  a  variable,  placed  in  an 
array,  or  written  to  a  permanent  data  file.  The  DMM  has  a 
continuous  data  sampling  rate  of  2.5  or  12.5  readings  per 
second.  This  sampling  rate  can  be  altered  as  desired  by  a 
command  from  the  BASICA  application  program.  The  slower 
sampling  rate  has  a  higher  accuracy  (+/-  .03%  of  the  input 
voltage)  than  the  higher  sampling  rate  (+/-  .05%  of  the 
input  voltage) .  Thus,  the  slower  sampling  rate  is  used 
throughout  the  experiments. 

When  a  data  reading  is  called  for,  the  last  complete 
reading  obtained  by  the  DMM  is  output  to  the  computer.  Speed 
of  data  acquisition  is  limited  to  the  speed  in  which  the  IBM 
computer  carries  out  each  BASICA  command.  This  varies 
between  15-120  milliseconds  per  command.  When  input  channels 
of  the  Relay  Multiplexer  are  switched  a  delay  of 
approximately  one  second  is  required  to  allow  the  circuit  to 
settle  before  a  measurement  is  taken.  The  result  is  a 
maximum  data  acquisition  rate  of  approximately  ten  readings 
per  second  when  one  channel  is  input  and  only  one  reading 
per  second  when  two  or  more  consecutive  channels  are  input. 
The  HP  data  acquisition  system  is  essentially  a  steady  state 
system. 
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3 .  Relay  Actuator 

The  Relay  Actuator  contains  eight  independently 
controllable  electrical  relays.  Each  of  the  eight  channels 
consists  of  one  single  pole-single  throw  relay.  Unlike  the 
Relay  Multiplexer,  all  eight  channels  can  be  set  in  any 
configuration  (open  or  closed)  simultaneously.  The  Relay 
Actuator  is  used  solely  in  controlling  the  angle  of 
incidence  setting  for  the  airfoil. 

4 .  Digital  Input  Port 

The  Digital  Input  Port  is  able  to  read  digitized 
data  in  four  different  modes:  Binary,  Decimal,  Hexadecimal, 
and  Octal. Each  of  these  modes  can  be  displayed  on  the 
computer  screen  when  the  Digital  Input  Port  is  interactive 
through  the  Soft  Front  Panels  software.  When  a  measurement 
is  called  from  a  BASICA  application  code  however,  the  output 
to  the  computer  is  always  in  the  Decimal  mode  with  2's 
complement  and  positive  logic.  In  this  mode  the  Digital 
Input  Port  assigns  a  bit  weight  to  each  of  the  input  pins 
that  read  "high"  from  the  Para-Scientific  Pressure  Computer. 
Digital  output  from  the  Para-Scientific  Pressure  Computer  is 
discussed  in  detail  in  Chapter  II.  D.  5. 

If  the  first  pin,  or  most  significant  bit,  reads 
high  the  Digital  Input  Port  assigns  a  negative  sign  to  the 
output.  A  typical  Decimal  value  corresponding  to  a  pressure 
value  of  14.73  PSI  looks  like  this:  -32766.  A  subroutine  in 
CALIB  and  PART2  inputs  the  Decimal  value  and  scales  it  to 
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determine  which  pins  read  high  and  which  give  no  input  (read 
low) .  When  scaling  is  complete  the  input  Decimal  value  is 
converted  to  its  numeric  equivalent  (in  this  example  -32766 
becomes  14.73)  in  pounds  per  square  inch. 

The  Digital  Input  Port  is  set  to  read  input 
digitized  data  continuously.  The  Para-Scientific  Pressure 
Computer  updates  its  pressure  reading  approximately  3  times 
per  second.  Therefore  a  new  pressure  reading,  in  digital 
form,  is  sent  to  the  Digital  Input  Port  3  times  per  second. 
When  a  measurement  is  called  from  the  BASICA  application 
code  of  CALIB  or  PART2 ,  the  last  complete  measurement  taken 
after  the  measure  statement  is  sent  to  the  computer.  This 
value  is  an  instantaneous  measurement  from  the  Para- 
Scientific  Transducer  (not  RMS) .  To  reduce  the  chance  of 
spurious  inputs,  70  measurements  are  called  for 
consecutively  and  averaged  in  the  computer  program.  Each 
measurement  takes  approximately  30  milliseconds  so  the  final 
pressure  reading  spans  6  pressure  updates  from  the  Para- 
Scientific  Pressure  Computer. 

5.   Para-Scientific  Pressure  Computer 

The  Para-Scientific  Pressure  Computer  can  accept 
inputs  from  up  to  eight  pressure  transducers.  Two  are 
utilized  for  data  acquisition  in  this  thesis  project.  The 
first  input  (PI)  is  from  a  differential  pressure  transducer. 
The  differential  transducer  reads  pitot-static  pressure  from 
the  individual  total  and  static  probes  in  the  test  section. 
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The  static  pressure  line  is  split  by  means  of  a  "T" 
connector.  One  line  is  connected  to  the  differential 
pressure  transducer,  the  other  is  connected  to  an  absolute 
pressure  transducer.  Static  pressure  from  the  absolute 
transducer  makes  up  the  second  input  (P2)  to  the  Para- 
Scientific  Pressure  Computer 

The  pressure  transducers  used  by  the  Pressure 
Computer  are  factory  calibrated.  Coefficients  are  supplied 
with  each  transducer  that  enable  the  Pressure  Computer  to 
convert  pressure  readings  to  an  output  of  PSI.  The  Pressure 
Computer's  output  is  updated  continuously  at  a  rate  of 
approximately  3  readings  per  second.  Output  data  can  be 
analog  or  digital  (BCD) .  Digital  output  is  used  in  this 
thesis  for  speed  and  accuracy. 

The  Pressure  Computer  provides  a  50  pin  digital 
output  to  accomodate  pressure  readings  as  high  as  10,000  PSI 
with  an  accuracy  to  the  fifth  decimal  (0.00001  PSI).  The 
limitation  on  this  system  was  in  the  HP  Digital  Input  Port. 
The  Digital  Input  Port  only  accepts  16  pins  of  input  data  so 
the  16  most  significant  data  pins  for  this  experiment  had  to 
be  used.  Data  taken  from  the  Pressure  Computer  and  input  to 
the  Digital  Input  Port  is  outlined  in  Figure  5.  The  accuracy 
of  this  configuration  is  +/-  0.002  PSI.  The  Pressure 
Computer  can  measure  pressures  from  its  two  input 
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Pin  Connections  to  the  HP  Digital  Input  Port 

Figure  5 

transducers  over  a  range  of  0.002  -  19.998  PSI.  The  scope  of 
possible  pressure  values  to  be  encountered  in  this  thesis 
project  are  well  within  this  measurement  range. 

E.   APPLICATION  SOFTWARE 

The  BASICA  application  program  that  runs  the  hotwire 
experiment  is  composed  of  three  parts.  The  first  is  called  a 
STATEFILE.  The  Statefile  is  a  program  automatically  compiled 
by  the  HP  Soft  Front  Panels  software.  It  tells  the  computer 
what  configuration  each  unit  was  left  in  when  the 
instruments  were  last  used  (i.e.,  Relay  Multiplexer  was  set 
for  channel  one  as  the  input  and  the  output  device  was 
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enabled,  etc.)-  The  Statefile  in  use  in  the  hotwire 
laboratory  is  named  HOTWIRE. HPC  and  is  located  in  the 
hotwire  sub-directory  on  the  computer's  fixed  C-disk. 

The  second  part  of  the  BASICA  application  program  is 
called  the  PROGRAM  SHELL.  The  Program  Shell  consists  of 
BASICA  program  lines  that  perform  initialization  chores  to 
allow  communication  between  the  HP  instruments  and  the  IBM 
computer.  In  essence,  the  Program  Shell  lets  the  computer 
know  every  function  available  at  each  of  the  data 
acquisition  instruments.  When  one  of  these  functions  is 
called  up  in  the  BASICA  application  program  the  computer 
already  "knows"  that  function  exists  and  where  to  find  it. 
The  Program  Shell,  like  the  Statefile,  is  created  by  the  HP 
system  Soft  Front  Panel  software. 

The  third  part  of  the  BASICA  application  program  is  the 
actual  application  code.  The  Program  Shell  occupies  lines  1 
through  999.  The  application  code  starts  at  line  1000  and 
begins  to  run  only  after  all  system  initialization  is 
complete.  The  application  code  consists  of  a  combination  of 
PC  instruments  statements  with  BASICA  keywords  and 
programming  statements.  The  hotwire  laboratory  is  run  in  two 
parts  using  BASICA  programs  "CALIB"  and  "PART2".  These 
programs  are  discussed  in  detail  in  Chapters  4  and  5, 
respectively . 

Certain  limitations  exist  when  using  the  BASICA 
language  to  run  a  long,  complicated  data  acquisition 
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program.  A  BASICA  application  program,  including  the  Program 
Shell  must  be  less  than  64  kilobytes  in  length.  This  is  a 
limitation  imposed  by  BASICA  version  2.0  in  use  by  the  IBM 
Disk  Operating  System.  Since  the  application  code  contains 
all  screen  graphics,  PC  instrument  statements,  plus  data 
manipulation  and  storage,  the  code  can  easily  exceed  this  64 
Kilobyte  limitation.  For  this  reason  the  hotwire  laboratory 
experiment  exists  in  two  parts. 

Transfer  of  control  from  CALIB  to  PART2  is  accomplished 
automatically  through  a  BASICA  "LOAD"  and  "RUN"  command. 
When  LOAD  PART2  is  executed  the  BASICA  program  CALIB  is 
removed  from  the  allotted  working  memory  and  PART2  is  put  in 
its  place.  The  RUN  command  automatically  causes  the 
experiment  program  "PART2"  to  run  without  any  further  user 
input.  When  employed  in  this  manner,  the  application  codes 
can  be  as  long  and  complicated  as  necessary  and  never  exceed 
the  memory  allocation  of  64  Kilobytes. 

Seven  named  variables  are  transferred  from  CALIB  to 
PART2.  Four  relate  the  calibrated  range  of  wind  tunnel 
velocities  for  subsequent  use  by  the  hotwire  experiment. 
Two  variables  contain  the  calibration  curve  slope  and 
intercept  data.  The  seventh  variable  contains  the 
calibration  tunnel  temperature  in  degrees  Farenheit.  This  is 
used  for  hotwire  compensation  when  the  tunnel  temperature 
varies  during  the  experiment.  The  seven  variables  are 
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written  to  a  permanent  data  file  that  is  updated  with  each 
experiment.  The  updated  data  file  is  read  into  PART2  after 
the  RUN  command  has  been  executed. 
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III.  DATA  ACQUISITION 

A.   ANGLE  OF  INCIDENCE 

A  ten  turn  potentiometer  is  located  on  the  drive 
mechanism  of  the  wind  tunnel  turntable.  Potentiometer  data 
is  routed  to  signal  amplifier  number  two  in  the 
instrumentation  rack.  Signal  amplifier  output  is  then 
directed  to  channel  five  of  the  Relay  Multiplexer.  Voltage 
output  from  the  signal  amplifier  is  zeroed  at  0  degrees  AOA 
and  is  positive  for  forward  turntable  motion  and  negative 
for  reverse  motion.  Voltage  readings  are  scaled  and 
converted  to  degrees  for  angle  of  incidence  measurements  in 
the  BASICA  application  code  of  CALIB  and  PART 2 . 

Turntable  control  for  forward  motion  (positive  AOA)  is 
routed  to  channel  one  of  the  Relay  Actuator.  Turntable 
control  for  reverse  motion  (negative  AOA)  is  routed  to 
channel  two  of  the  Relay  Actuator.  Figure  6  shows  signal 
path  channels  leading  to  turntable  control.  When  forward 
turntable  motion  is  desired,  Relay  Actuator  channel  one 
contacts  are  closed.  Timing  for  Relay  Actuator  closure  and 
turntable  motion  is  scaled  to  account  for  the  diminishing 
difference  between  the  AOA  desired  and  the  present  AOA.  When 
the  turntable  nears  the  desired  AOA  the  relay  contacts  will 
close  for  only  a  portion  of  a  second  to  allow  for  an 
accuracy  of  +/-  0.1  degree.  After  each  Relay  Actuator 
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closure  (resulting  in  turntable  motion)  the  present  AOA  is 
read  and  the  new  difference  between  the  desired  AOA  and 
present  AOA  is  computed.  This  sequence  is  identical  for 
reverse  turntable  motion  except  Relay  Actuator  channel  two 
contacts  close. 

B.  TRAVERSE  ASSEMBLY  POSITION 

A  potentiometer  on  the  traverse  assembly  sends  its 
output  to  signal  conditioner  number  three,  located  in  the 
top  of  the  instrumentation  rack.  The  potentiometer  rides  a 
gear  train  across  the  wind  tunnel  test  section. 
Potentiometer  output  is  zeroed  at  the  near  end  of  the  test 
section  and  increases  linearly  to  the  far  side.  A  DC  voltage 
signal  from  signal  conditioner  number  three  is  routed  to 
input  channel  four  of  the  Relay  Multiplexer.  Scaling  of  the 
input  voltage  to  "inches"  takes  place  in  the  BASICA 
application  code  of  PART2.  Figure  7  shows  the  traverse 
assembly  mounted  in  the  wind  tunnel  test  section  behind  the 
airfoil  model. 

C.  PITOT-STATIC  PRESSURE 

The  upper  sensor  on  the  traverse  assembly  in  Figure  8 
senses  wind  tunnel  total  pressure.  This  pressure  makes  up 
one  half  of  a  pitot-static  reading  for  tunnel  flow  field 
measurement  and  is  connected  to  one  side  of  a  differential 
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pressure  transducer.  The  longest  sensor  on  the  traverse  unit 
(third  from  the  top  in  Figure  8)  senses  static  pressure. 
This  reading  makes  up  the  second  half  of  the  pitot-static 
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Pitot-Static/Hotwire  Traverse  assembly 

Figure  8 

measurement  and  is  routed  to  the  other  side  of  the 
differential  pressure  transducer.  Pressure  data  from  the 
differential  transducer  is  routed  to  input  .PI  on  the  Para- 
Scientific  Pressure  Computer. 

Although  the  two  sensors  occupy  separate  tubes,  they 
function  exactly  as  a  conventional  pitot  tube.  Signal  paths 
for  this  system  are  shown  in  Figure  9 .  The  traverse  assembly 
pitot-static  system  was  constructed  in  this  manner  to 
measure  static  and  stagnation  pressures  at  the  same  axial 
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location.  The  total  and  static  pressure  sensors  as  well  as 
the  hotwire  probe,  all  obtain  their  flow  measurements  on  the 
same  vertical  plane  aft  of  the  airfoil.  Wake  profile 
readings  taken  from  the  hotwire  or  pitot-static  system  are 
virtually  identical. 

Static  pressure  in  the  test  section  is  also  used  for 
air  density  and  flow  velocity  calculations.  The  static 
pressure  line  that  is  routed  to  the  differential  pressure 
transducer  is  also  connected  to  an  absolute  pressure 
transducer  by  means  of  a  "T"  connection.  Static  pressure 
readings  from  the  absolute  pressure  transducer  are  connected 
to  input  P2  on  the  Para-Scientific  Pressure  Computer. 

D.  HOTWIRE  VOLTAGE  DATA 

Hotwire  voltage  data  is  output  directly  in  DC  volts  from 
the  back  of  the  IFA  100  thermal  anemometer.  The  voltage 
output  is  routed  to  input  channel  1  of  the  Relay 
Multiplexer.  This  value  is  read  directly  by  the  Digital 
Multi-Meter  and  at  no  time  are  the  values  scaled  or  altered. 

E.  WIND  TUNNEL  TEMPERATURE 

A  Copper-Constantin  thermocouple  is  located  inside  the 
wind  tunnel  in  the  converging  section  prior  to  the  test 
section.  Thermocouple  output  is  routed  directly  to  input 
channel  seven  of  the  Relay  Multiplexer.  An  internal  sensor 
in  the  Relay  Multiplexer  monitors  and  converts  ambient  room 
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temperature  to  a  reference  output  voltage.  This  reference 
voltage  is  connected  to  input  channel  eight  of  the  Relay- 
Multiplexer.  Voltage  output  from  the  temperature 
thermocouple  to  channel  7  and  reference  voltage  output  to 
channel  8  are  read  by  the  DMM  consecutively  and  their  values 
are  added  to  obtain  a  final  accurate  wind  tunnel  temperature 
reading.  A  subroutine  for  measuring  wind  tunnel  temperature 
is  located  in  both  CALIB  and  PART2 . 

NOTE:  Coefficient  values  used  in  the  temperature  measurement 
subroutine  are  extremely  dependent  on  the  type  of 
temperature  thermocouple  used.  If  the  present  Copper- 
Constantin  thermocouples  are  replaced  or  altered  in  future 
research  the  HP  PC  Instrument  User's  Manual  should  be 
consulted  for  information  on  new  thermocouple  coefficients. 
Any  change  to  the  present  system  can  result  in  spurious  or 
inaccurate  tunnel  temperature  readings. 

NOTE:  Figure  10  contains  all  signal  paths  utilizing  the  HP 
Relay  Multiplexer.  Data  acquisition  channels  not  using  the 
HP  Relay  Multiplexer  were  illustrated  in  Figures  6  and  9. 
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IV.  HOTWIRE  CALIBRATION 

Heat  loss  from  the  hotwire  sensor  due  to  forced 
convection  is  dependent  on  many  parameters: 

flow  velocity  -  magnitude  and  direction 
probe  geometry  and  material 
temperature  of  the  sensor 
Although  some  of  these  influences  can  be  kept  constant  for  a 
particular  experiment  no  sufficiently  exact  theoretical 
relationship  for  thermal  anemometer  heat  transfer  exists. 
Thus,  the  need  for  system  calibration. 

For  a  given  probe  in  the  low  speed  wind  tunnel  the  wire 
response  E,  for  constant  temperature  anemometry  is  a 
function  of: 

flow  velocity 
yaw  angle 
pitch  angle 
or  E  =  f (velocity,  yaw  angle,  pitch  angle) .  The  yaw  and 
pitch  angle  remain  zero  in  this  thesis  as  the  traverse 
assembly  is  not  built  to  swivel.  The  probe  always  faces 
directly  into  the  flow  as  the  traverse  unit  moves  across  the 
test  section.  The  aim  of  hotwire  system  calibration  is  to 
determine  the  function  f  in  analytic  or  tabular  form. 

Hotwire  system  calibration  is  accomplished  through  the 
first  BASICA  program  called  CALIB.  To  calibrate  a  hotwire 
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sensor,  a  comparison  is  made  of  the  hotwire  voltage  output 
to  a  known  value  of  flow  velocity  in  the  wind  tunnel  test 
section.  Air  flow  in  the  NPS  wind  tunnel  can  be  considered 
low  speed,  steady,  incompressible  flow.  Under  these 
conditions  a  pitot  tube  will  give  extremely  accurate  flow 
velocity  measurements.  Thus,  the  pitot  tube  becomes  the 
standard  by  which  hotwire  probe  output  is  compared. 

The  calibration  curve  is  created  by  running  the  wind 
tunnel  at  the  lowest  velocity  to  be  used  in  the  experiment 
and  incrementally  increasing  the  tunnel  velocity.  At  each 
discrete  wind  speed,  hotwire  output  and  tunnel  velocity 
values  are  sampled.  The  selected  tunnel  speeds  range  from 
the  lowest  to  highest  velocity  to  be  encountered  in  the 
experiment.  Since  the  hotwire  system's  accuracy  is  only 
known  within  the  calibration  range-,  all  expected 
experimental  velocities  are  covered.  CALIB  uses  a  least 
squares  analysis  to  determine  the  best  linear  fit  through 
the  calibration  data  points.  Figure  11  is  representative  of 
all  calibration  curves  using  this  method.  Deviation  from  the 
linear  fit  for  each  data  point  is  typically  very  small. 

One  suggestion  for  future  experiments  at  NPS  is  to  run 
the  wind  tunnel  during  the  data  collection  phase  of  an 
experiment  in  the  middle  of  the  calibrated  range  of 
velocities.  If  a  series  of  velocities  are  to  be  used,  care 
should  be  taken  to  ensure  that  they  are  within  the 
calibration  range.  A  calibration  curve  can  be  constructed  at 


36 


LU 
> 

co 


o 


o 


< 

cr 

CD 


< 

o 

LU 

cr 


o 

IE 


o 
in 


o 
cq 

IT) 


o 
in 

in 


o 

in 


o 
ro 

in 


o 

in 


o 
.to 

IOlO 

d 

UJ 

°> 

o 


I  M  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  !  I  I  I  M  I  I  I  I  I  I  I  I  I  I  I  I  I  I  U  I 


O 


oo  LU 
600 


o 

o 


.SLiJ 


o 


> 


o 

°  uZ 

ciO 
'-LU 

t    cyo 


■Sen 

dLU 


f-O 

00 


O 

in 


CD 


(SilOA)    0"23    30V110A   BdlMlOH 


3 

oc 


37 


any  time.  If  the  validity  of  experimental  data  is  in 
question  it  is  best  to  go  back  to  the  begining  and  re-create 
the  calibration  curve. 

A  minimum  of  10  data  points  should  be  taken  to 
construct  the  calibration  curve.  The  accuracy  of  the  curve 
is  directly  related  to  the  number  of  inputs.  For  time 
consideration  in  classroom  demonstrations,  10  points  will 
prepare  a  curve  with  accuracy  to  1-2%.  If  the  hotwire  system 
is  to  be  used  for  research  or  other  important  measurements, 
20-3  0  data  points  are  necessary  to  achieve  optimum  accuracy. 
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V.  DATA  ANALYSIS 

A.   EXPERIMENTAL  TECHNIQUE 

The  experimental  portion  of  this  thesis  is  performed  by 
the  BASICA  application  code  called  PART 2 .  PART2  starts  by- 
prompting  the  user  to  input  desired  angles  of  attack  to  be 
analyzed.  When  the  first  angle  of  attack  is  set  in  the  test 
section  a  flow  field  velocity  update  is  performed.  If 
adjustments  are  necessary  the  user  can  re-adjust  test 
section  wind  speed  and  update  the  tunnel  velocity 
repeatedly. 

After  the  desired  angle  of  incidence  and  test  section 
flow  velocity  are  set,  the  system  is  ready  for  data 
sampling.  Traverse  position  in  the  test  section  is  initially 
displayed  on  the  computer  screen.  The  researcher  can  move 
the  traverse  assembly  anywhere  across  the  test  section  for 
the  first  data  point.  Accuracy  of  the  traverse  position  is 
+/-  0.01  inch. 

The  airfoil  wake  velocity  profile  is  the  focus  of 
interest  in  this  experiment.  In  order  to  concentrate  on  the 
wake  profile  and  not  spend  excessive  time  sampling  the  free 
stream,  correlation  with  the  oscilloscope  is  necessary.  When 
the  hotwire  system  is  running,  its  electrical  response  to 
the  test  section  flow  field  is  output  to  the  oscilloscope. 
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Oscilloscope  Response  To  Steady  Flow 
Figure  12 

The  oscilloscope  trace  will  be  stable  under  steady  flow 
conditions  (see  Figure  12) .  As  the  hotwire  sensor  nears  the 
airfoil  wake  the  oscilloscope  trace  will  begin  to  jump  in 
response  to  the  disturbed  air  flow  at  the  wake  boundary  (see 
Figure  13) .  When  the  hotwire  sensor  is  in  the  fully 
developed  wake  the  oscilloscope  will  display  the  unsteady 
response  caused  by  turbulent  flow  (see  Figure  14) .  By 
correlating  the  hotwire  sensor  response  as  seen  on  the 
oscilloscope  with  the  traverse  position  as  given  on  the 
computer  screen,  the  start  and  end  of  the  wake  boundary  can 
easily  be  detected. 
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Oscilloscope  Response  To  Wake  Boundary 

Figure  13 


Oscilloscope  Response  To  Turbulent  Flow 

Figure  14 
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Experimental  procedure  consisted  of  moving  the  traverse 
assembly  toward  the  wake  until  an  unsteady  response  was 
displayed  on  the  oscilloscope.  Once  the  position  of  the  wake 
edge  is  known,  the  traverse  unit  is  backed  out  of  the  wake 
into  the  freestream  and  measurements  are  initiated.   By 
correlating  traverse  unit  position  with  hotwire  sensor 
response,  experimental  data  sampling  can  be  accomplished  in 
a  more  efficient  manner. 

When  the  traverse  assembly  is  positioned  for  a  reading, 
the  computer  samples  the  hotwire  system  voltage  and  pitot- 
static  pressures  for  independent  test  section  velocity 
measurements.  This  procedure  is  repeated  as  many  times  as 
desired  to  obtain  a  full  wake  velocity  profile  of  the 
airfoil.  If  PART2  is  to  be  used  as  a  classroom  demonstration 
as  few  as  three  data  points  per  inch  of  traverse  travel 
across  the  test  section  will  provide  a  viable  velocity 
profile  graphical  output.  If  the  experiment  is  to  include  a 
drag  coefficient  measurement  of  the  airfoil,  at  least  five 
data  points  per  inch  of  traverse  travel  is  recommended.  Ten 
points  per  inch  is  optimal. 

B.   DRAG  COEFFICIENT  MEASUREMENT 

PART2  is  capable  of  measuring  the  drag  coefficient  of 
the  airfoil  at  each  angle  of  incidence.  The  drag  coefficient 
is  deduced  as  a  loss  of  linear  momentum  in  the  air  flow  as 
it  passes  over  the  airfoil.  At  each  data  sampling  position 
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the  hotwire  and  pitot-static  readings  are  converted  into  a 
dynamic  pressure  value.  Dynamic  pressure  in  the  wake,  q,  is 
compared  in  a  ratio  to  the  dynamic  pressure  of  the 
freestream,  qo.  The  dynamic  pressure  ratio  q/qo,  is 
approximated  by  the  following  equation: 

q/qo  =  V2/V1(1-V2/V1) 
where  V2  is  the  flow  velocity  in  the  wake  and  VI  is  the 
freestream  velocity. 

This  method  of  analyzing  the  wake  velocity  profile  will 
only  work  if  the  traverse  unit  is  situated  far  enough  behind 
the  wing  so  that  the  wake  has  returned  to  tunnel  static 
pressure.  In  early  experimental  trials  of  this  thesis  the 
static  pressure  profile  was  measured  across  the  wake  for  all 
angles  of  incidence.  In  every  case  the  change  in  static 
pressure  across  the  test  section  was  zero.  Thus,  obtaining 
profile  drag  coefficients  by  the  momentum  method  in  the  NPS 
wind  tunnel  is  fully  justified  (see  Chapter  VI) .  Stalled 
conditions  were  not  included  in  this  analysis  as  the 
momentum  method  breaks  down.  A  large  part  of  the  drag  in 
this  case  is  caused  by  rotational  losses  and  does  not  appear 
as  a  drop  in  linear  momentum. 

The  dynamic  pressure  ratio,  q/qo,  is  calculated  for 
each  data  point  in  the  wake  velocity  profile.  If  the  dynamic 
pressure  ratio  were  to  be  plotted  on  an  X-Y  axis  (test 
section  position  on  the  X  axis,  q/qo  on  the  Y  axis)  it  would 
look  like  an  inverted  velocity  profile.  The  area  under  this 
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"curve",  multiplied  by  (2/airfoil  chord)  makes  up  an 
uncorrected  value  for  the  profile  drag  coefficient.  The  term 
uncorrected  refers  to  the  fact  that  the  drag  coefficient  is 
based  on  a  dynamic  pressure,  q,  which  would  exist  if  no 
airfoil  were  present  in  the  test  section.  The  uncorrected 
drag  coefficient  is  calculated  using  hotwire  system  data  as 
well  as  pitot-static  system  data  for  comparison  purposes. 

Integration  of  the  dynamic  pressure  ratio  is 
accomplished  by  a  modified  Trapezoidal  routine.  [Ref.  1: 
pp.  316-325]  Normally  trapezoidal  integration  is  carried  out 
by  averaging  a  function  over  an  incremental  area,  then 
multiplying  the  average  by  the  interval  length.  Round  off 
errors  can  become  excessive  if  the  function  to  be  integrated 
varies  rapidly.  To  reduce  round  off  errors  an  Akima  Spline 
fitting  technique  is  used  to  fit  a  third  order  polynomial 
between  each  pair  of  data  points.  Trapezoidal  integration  is 
then  used  to  sum  up  all  incremental  areas  with  a  cubic  fit 
between  data  pairs  vice  an  average  value. 

C.   TUNNEL  CORRECTIONS  FOR  DRAG  ANALYSIS 

Testing  an  airfoil  section  in  a  wind  tunnel  is  not  the 
same  as  testing  in  a  free  air  stream.  The  test  section  walls 
provide  constraints  to  the  air  flow  that  must  be  corrected. 
The  NACA  66 (215) -216  airfoil  tested  in  this  thesis  is  a  two 
dimensional  airfoil  section,  therefore  two  dimensional 
corrections  apply.  [Ref.  2: pp.  7.1-8.3] 
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dimensional  airfoil  section,  therefore  two  dimensional 
corrections  apply.  [Ref.  2: pp.  7.1-8.3] 

For  the  testing  conditions  encountered  in  this 
experiment  only  two  corrections  are  relevant.  The  first  is 
known  as  "solid  blocking".  With  the  airfoil  in  the  wind 
tunnel  the  test  section  cross  sectional  area  is  reduced.  A 
reduction  in  area  causes  the  surrounding  flow  to  speed  up 
according  to  Bernoulli's  law.  Solid  blocking  results  in  an 
increase  in  forces  and  moments  on  the  airfoil  at  any  given 
angle  of  incidence. 

The  second  correction  is  called  "wake  blocking".  The 
wake  behind  the  airfoil  has  a  lower  velocity  than  that  of 
the  freestream.  To  satisfy  continuity  the  air  flow  outside 
the  wake  must  speed  up  to  maintain  a  constant  volume  of  air 
flowing  through  the  test  section.  The  increase  in  free 
stream  velocity  near  the  wake  also  results  in  an  increase  in 
forces  and  moments  on  the  airfoil.  In  this  case  the  effect 
increases  with  an  increase  in  wake  size  (i.e.,  at  high 
angles  of  incidence) . 

The  solid  blocking  and  wake  blocking  corrections  are 
applied  to  the  uncorrected  drag  coefficients  to  produce  a 
final  corrected  profile  drag  coefficient. 
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VI.  RESULTS 

A.   WAKE  VELOCITY  PROFILES 

Figures  15  through  19a  represent  wake  velocity  profiles 
obtained  through  the  hotwire  and  pitot-static  systems.  In 
all  cases  the  profiles  are  remarkably  similar.  The  profiles 
highlight  two  main  points  concerning  hotwire  anemometry  in 
the  NPS  wind  tunnel.  First,  the  hotwire  anemometry  system 
works  quite  well.  Second,  when  turbulence  intensity  is  not 
too  high  (angles  of  incidence  less  than  ten  degrees)  the 
hotwire  system  has  the  greater  sensitivity  in  unsteady  flow. 
In  Figures  15,16,  and  17  the  hotwire  profile  consistently 
reaches  a  lower  peak  at  the  wake  center  than  the  pitot- 
static  system.  The  increased  sensitivity  falls  off  at  higher 
turbulence  intensities  caused  by  large  angles  of  incidence 
(see  Figures  19  and  19a) . 

An  explanation  of  the  spurious  data  points  or  jumps  in 
an  otherwise  smooth  velocity  profile  requires  a  reiteration 
of  how  the  data  is  obtained.  The  Digital  Multi-Meter  is 
capable  of  updating  its  voltage  measurements  (with  an 
acceptable  accuracy)  at  a  rate  of  only  2.5  readings  per 
second.  When  a  measurement  is  called  for  from  the 
application  programs,  the  last  full  measurement  obtained  by 
the  Digital  Multi-Meter  is  sent  to  the  computer.  These  "snap 
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shot"  measurements  reflect  the  unsteadiness  of  the  flow 
field  in  the  wake  of  the  airfoil. 

The  data  acquisition  system  used  in  this  thesis  was 
designed  to  measure  steady  flow  conditions.  To  simulate  a 
steady  state  (or  time  averaged)  data  reading  many 
consecutive  "snap  shot"  measurements  are  taken  by  the 
Digital  Multi-Meter  and  averaged  to  represent  a  mean  value. 
In  performing  a  wake  profile  analysis  a  trade  off  must  be 
accepted  between  the  number  of  measurements  taken  for  each 
data  point  and  the  time  allotted  for  data  collection. 

Each  of  the  velocity  profiles  except  Figure  19a  were 
prepared  using  70  consecutive  measurements  per  data  point. 
In  a  high  turbulence  intensity  flow,  70  measurements  per 
data  point  is  clearly  insufficient  for  a  smooth  graphical 
velocity  profile  output.  Figure  19  demonstrates  that  only  a 
rough  approximation  of  the  actual  velocity  profile  can  be 
made  in  high  turbulence  regimes.  The  hotwire  velocity 
profile  of  Figure  19a  was  prepared  using  300  measurements 
per  data  point  for  comparitive  purposes.  The  pitot-static 
profile  remained  at  70  readings  per  data  point.  The  hotwire 
curve  is  much  smoother  in  this  case,  proving  that  a  time 
averaged  value  of  the  wake  velocity  will  provide  the 
smoothest  plot.  It  is  also  interesting  to  note  that  the  wake 
profiles  measured  in  Figure  19  were  completed  in 
approximately  3  5  minutes.  The  wake  profiles  prepared  for 
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Figure  19a  required  over  90  minutes  for  completion  due  to 
the  extra  measurements  in  the  hotwire  velocity  profile. 

The  data  acquisition  system  utilizing  the  HP  Digital 
Multi-Meter  has  proven  to  be  a  useful  and  effective  system 
in  steady  flow  and  low  turbulence  intensity  regimes.  Outside 
of  these  flow  regimes  the  system  loses  its  effectiveness  and 
performance  becomes  marginal.  The  user  should  be  well  aware 
of  the  limitations  on  the  system.  Until  another  voltage 
measuring  device  capable  of  time  averaging  can  be 
incorporated  into  the  data  acquisition  system,  wake  profiles 
of  10  degrees  angle  of  incidence  and  above  should  be  viewed 
as  little  more  than  approximations.  Highly  accurate  and 
repeatable  profiles  can  still  be  obtained  in  low  turbulence 
intensity  flow;  i.e.,  angles  of  incidence  between  -9  and  +9 
degrees . 

B.   DRAG  COEFFICIENT  ANALYSIS 

Figure  20  displays  corrected  drag  coefficients 
calculated  from  the  hotwire  system  data  and  pitot-static 
system  data.  Both  plots  follow  the  same  general  shape, 
exhibiting  drag  buckets  and  steep  rises  in  Cd  at  high  lift 
coefficients.  It  should  be  noted  that  the  drag  coefficients 
for  both  systems  are  closer  in  value  at  lower  lift 
coefficients.  This  is  in  keeping  with  the  discussion  on 
system  limitations  in  part  A.  The  wake  velocity  profiles  are 
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are  generally  more  accurate  at  lower  angles  of  incidence, 
creating  more  accurate  and  repeatable  drag  coefficients. 

The  NACA  drag  coefficient  data  was  obtained  using  the 
same  wake  survey  method  as  used  in  this  thesis.  The 
difference  between  NPS  wind  tunnel  derived  drag  coefficients 
and  the  NACA  data  can  be  explained  in  several  ways.  The  NACA 
data  on  the  66 (215) -216  airfoil  was  collected  at  three 
discrete  Reynolds  number  flow  velocities,  the  slowest  being 
Re  =  3.0xE6.  Data  was  collected  in  the  NPS  wind  tunnel  at  a 
Reynolds  number  of  approximately  Re  =  6.0xE5.  The  slow 
tunnel  velocity  was  primarily  chosen  to  preserve  the  hotwire 
sensor  for  multiple  data  collection. 

The  airfoil  model  used  to  provide  the  NACA  data  had  a  2 
foot  chord  and  spanned  a  3  foot  wide  test  section.  The 
airfoil  model  used  in  the  NPS  wind  tunnel  has  a  1  foot  chord 
and  is  32  inches  in  length.  NACA  is  careful  to  point  out 
that  changes  in  surface  roughness,  Reynolds  number,  and 
model  geometry  will  affect  final  drag  coefficient 
calculations.  Therefore  an  exact  comparison  of  the  NACA  drag 
coefficient  data  with  NPS  wind  tunnel  derived  data  is 
unproductive.  A  comparison  of  trends  is  better  suited  for 
the  two  differing  test  facilities.  One  can  conclude  that  the 
wake  survey  method  used  in  this  thesis  will  provide  a  very 
accurate  value  of  the  drag  coefficient  for  the  NPS  66(215)- 
216  model  airfoil. 
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VII.   CONCLUDING  REMARKS 

A  data  acquisition  system  has  been  installed  in  the  NPS 
wind  tunnel  that  utilizes  hotwire  anemometry  to  analyze  the 
wake  velocity  profiles  of  two  dimensional  airfoil  models.  A 
conventional  pitot-static  system  is  also  incorporated  in  the 
data  acquisition  process  for  comparative  purposes.  The 
hotwire  data  acquisition  system  is  self  contained  and  needs 
no  interface  with  any  other  experiment  currently  in  use  in 
the  NPS  wind  tunnel. 

The  hotwire  data  acquisition  system  is  designed  in  two 
separate  parts  to  enhance  versatility.  System  calibration 
can  be  performed  as  a  classroom  demonstration  or  as  the 
prelude  to  an  entire  wake  profile  analysis.  The  wake  profile 
analysis  can  also  be  used  as  a  classroom  demonstration  or  as 
a  complete  experimental  laboratory.  In  either  case,  all 
graphical  and  tabular  data  are  available  for  online  computer 
display  or  hardcopy  output. 

A  NACA  66 (215) -216  airfoil  model  was  analyzed  to  verify 
system  performance.  Virtually  any  two  dimensional  wing 
section  that  will  fit  in  the  wind  tunnel  test  section  can  be 
used  for  research.  Wake  velocity  profiles  and  drag 
coefficients  can  also  be  found  for  non-airfoil  shaped 
models.  The  lessons  learned  from  this  thesis  project  will  be 
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of  value  to  aeronautical  laboratory's  held  at  NPS  as  well  as 
future  research  and  improvements  to  the  Low  Speed  Wind 
Tunnel  Facility. 
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VIII.   RECOMMENDATIONS 

This  thesis  project  was  designed  and  built  as  a  steady 
state  data  acquisition  system.  The  next  logical  step  for 
hotwire  anemometry  research  is  to  take  advantage  of  its  fast 
response  time.  Hardware  and  software  additions  should  be 
made  so  the  IBM  PC-AT  computer  can  handle  data  inputs  of  up 
to  100  kilohertz.  With  these  changes  installed,  a  real  time 
analysis  of  wake  turbulence,  velocity  vector  components,  and 
turbulence  intensity  can  be  realized.  One  short  term 
improvement  that  would  greatly  enhance  current  measurement 
capabilities  is  to  replace  the  Hewlett  Packard  Digital 
Multi-Meter  with  a  voltmeter  capable  of  outputting  a  mean 
voltage  reading.  With  this  change  incorporated,  all  velocity 
profiles  up  to  the  stall  region  would  be  smooth  and 
continuous. 

Hewlett  Packard  now  produces  a  software  package  that 
allows  their  data  acquisition  equipment  to  be  compatible 
with  Quick  Basic  (a  compiled  form  of  Basic) .  Quick  Basic  is 
approximately  ten  times  faster  in  execution  than  the 
Advanced  Basic  (BASICA)  used  in  this  thesis  project.  If 
CALIB  and  PART2  were  executed  in  Quick  Basic  then  hundreds 
of  measurements  per  data  point  would  become  a  reality  since 
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time  constraints  would  become  negligible.  This  would 
eliminate  the  need  for  a  voltmeter  with  mean  value  output 
capability. 

The  traverse  assembly  is  one  of  the  greatest 
limitations  to  experimental  efficiency.  The  current 
configuration  requires  two  people  to  coordinate  traverse 
movement  through  the  use  of  hand  signals  (to  compensate  for 
tunnel  noise)  or  one  person  to  repeatedly  move  back  and 
forth  between  the  wind  tunnel  test  section  and  the  computer. 
A  motorized  traverse  assembly  that  can  be  computer 
controlled  with  no  loss  in  accuracy  is  highly  recommended. 

Another  system  deficiency  resides  in  the  turntable 
positioning  unit  for  angle  of  incidence  control.  The  chain 
drive  assembly  has  play  in  its  movement  that  cannot  be 
compensated  for  in  a  computer  program.  The  relay/ DC  motor 
system  is  imprecise  and  must  rely  on  software  timing  for 
correct  positioning.  The  signal  amplifier  will  waiver  as 
much  as  0.2  degrees  after  turntable  motion  has  ceased.  To 
compensate  for  this  lack  of  precision  it  is  recommended  that 
angles  of  incidence  be  entered  into  PART2  one  at  a  time  if 
multiple  experiments  are  to  be  performed.  As  long  as  the 
tunnel  or  hotwire  system  is  not  shut  down,  PART2  can  be 
restarted  repeatedly  without  altering  system  calibration. 

The  Hewlett  Packard  Laser  Jet  Printer  also  contributes 
to  inefficiency  in  the  laboratory.  Graphical  outputs  each 
take  between  4-6  minutes  for  hardcopy  output.  During  this 
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time  the  computer  is  continually  sending  data  to  the  printer 
and  is  not  available  for  other  use.  A  faster  printer  or  a 
printer  buffer  that  would  free  the  computer  for  other  uses 
would  eliminate  unproductive  waiting  during  classroom 
demonstrations.  This  upgrade  is  desireable  for  laboratory 
work  as  well  and  is  highly  recommended. 

System  acquisition  for  this  thesis  project  seems  to 
have  been  performed  in  a  sporadic  manner.  Large  amounts  of 
money  were  invested  in  top  quality  equipment  (ie,  the  HP 
Data  Acquisition  System  and  the  Para-Scientific  Pressure 
Computer)  but  their  performance  becomes  marginal  when  used 
as  a  cohesive  data  acquisition  unit.  For  example,  the  Para- 
Scientific  Pressure  Computer  provides  a  50  pin  digital 
output  to  supply  pressure  data  with  a  +/-  0.00001  PSI 
accuracy.  The  HP  Digital  Input  Port  can  only  accept  16  data 
pins,  therefore  the  precision  of  the  Pressure  Computer  is 
now  reduced  to  +/"  0.002  PSI.  This  accuracy  is  only  a 
marginal  improvement  over  the  analog  pressure  transducers 
previously  installed  in  the  wind  tunnel  test  section.  It  is 
highly  recommended  that  any  new  equipment  purchased  for 
upgrading  present  wind  tunnel  capabilities  be  carefully 
scrutinized  for  compatibility  with  equipment  already  in  use. 
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APPENDIX  A 

Sample  Calculations  of  Test  Section  Correction  Factors 

Solid  Blocking.  Esb 

Esb  =  Kl  *  Model  Volume 
CA1.5 

where  Kl  =  0.52  for  a  model  spanning  the  test  section 

vertically 

C  =  tunnel  test  section  area  (8.75  ftA2) 

Model  Volume  =  .7  *  max  thickness  *  chord  *  span 

Model  Volume  =  451.584  inA3 

Model  Volume  =  .2613  ftA3 


Esb  =  .52  *  .2613 
8,75A1.5 

Esb  =  0.00525 


Wake  Blocking.  Ewb 
Ewb  =  (c/2W)  *  Cdu 
where  c  =  model  chord  =  12  in 

W  =  test  section  width  =  45  in 

Cdu  =  uncorrected  drag  coefficient  determined  for 
each  experimental  profile 

Corrected  Drag  Coefficient.  Cd 
Cd  =  Cdu  *  (1-  3Esb  -  2Ewb) 
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APPENDIX  B 


Program  Listing:  Shell  Program  and  CALIB 


1  DEF  SEG: CLEAR  ,&HFE00:  KEY  OFF:  COLOR  1,7,7:  LOCATE  12,26: 
PRINT  "INITIALIZING  HP  INSTRUMENTS":  GOTO  4 

2  GOTO  1000  '  User  program 

3  GOTO  900  '  Error  handling 

4  I=&HFE00  '  Copyright  Hewlett-Packard  1984,1985 

5  PCIB.DIR$=ENVIRON$("PCIB") :  CLS :  COLOR  1,7,7:  LOCATE 
12,26:  PRINT  "INITIALIZING  HP  INSTRUMENTS" 

6  I$=PCIB.DIR$+"\PCIBILC.BLD" 

7  BLOAD  I$,I 

8  CALL  I(PCIB.DIR$,I%,J%) :PCIB.SEG=I% 

9  IF  J%=0  THEN  GOTO  13  * 

10  PRINT  "Unable  to  load."; 

11  PRINT  "    (Error  #";J%;")" 

12  END 

13  ' 

14  DEF  SEG=PCIB.SEG:O.S=5:C.S=10:I.V=15 

15  I.C=20:L.P=25:LD.FILE=30 

16  GET.MEM=35:L.S=40:PANELS=4  5:DEF.ERR=50 

17  PCIB.ERR$=STRING$(64,32)  :  PCIB.NAME$=STRING$ (16 , 32) 

18  CALL  DEF.ERR(PCIB.ERR,PCIB.ERR$,PCIB.NAME$,PCIB.GLBERR)  : 
PCIB.BASERR=255 

19  ON  ERROR  GOTO  3 
2  0  J=-l 

21  I$=PCIB.DIR$+"\PCIB.SYN" 

22  CALL  O.S(I$) 

23  IF    PCIB.ERRO0    THEN    ERROR   PCIB.BASERR 

24  1=0 

25  CALL 

I . V ( I , READ . REGISTER , READ . SELFID , DEFINE , INITIALIZE . SYSTEM) 

2  6  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

27  CALL 

I . V ( I , ENABLE . SYSTEM , DISABLE . SYSTEM , INITIALI ZE , POWER . ON ) 

28  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

29  CALL  I. V( I, MEASURE, OUTPUT, START, HALT) 

30  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

31  CALL 

I . V ( I , ENABLE . INT . TRIGGER , DISABLE . INT . TRIGGER , ENABLE . OUTPUT , D 
I SABLE. OUTPUT) 

32  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

3  3  CALL  I. V( I, CHECK. DONE, GET. STATUS, SET. FUNCTION, SET. RANGE) 
34  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

3  5  CALL  I. V( I, SET. MODE, WRITE. CAL, READ. CAL, STORE. CAL) 
3  6  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

37  CALL  I. V( I, DELAY, SAVE. SYSTEM, J, J) 

38  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 
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39  1=1 

40  CALL  I. V( I, SET. GATETIME, SET. SAMPLES, SET. SLOPE, SET. SOURCE) 

41  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

42  CALL  I. C( I, FREQUENCY, AUTO. FREQ, PERIOD, AUTO. PER) 

43  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

44  CALL  I. C( I, INTERVAL, RATIO, TOTALIZE,R100MILLI) 

45  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

46  CALL  I.C(I,R1,R10,R100,R1KILO) 

47  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

48  CALL  I.C(I,R10MEGA,R100MEGA,CHAN.A,CHAN.B) 

49  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

50  CALL  I. C( I, POSITIVE, NEGATIVE, COMN, SEPARATE) 

51  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

52  1=2 

53  1=3 

54  CALL  I. V(I, ZERO. OHMS, SET. SPEED, J, J) 

55  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

56  CALL  I.C(I,DCVOLTS,ACVOLTS,OHMS,R2  00MILLI) 

57  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

58  CALL  I.C(I,R2,R20,R200,R2KILO) 

59  IF   PCIB.ERRO0    THEN   ERROR   PCIB.BASERR 

60  CALL   I.C(I,R20KILO,R200KILO,R2MEGA,R2  0MEGA) 

61  IF   PCIB.ERRO0    THEN   ERROR   PCIB.BASERR 

62  CALL   I.C(I,AUTOM,R2.5,R12.5,J) 

63  IF    PCIB.ERRO0    THEN   ERROR   PCIB.BASERR 

64  1=4 

65  CALL 

I . V ( I , SET . COMPLEMENT , SET . DRIVER , OUTPUT . NO . WAIT , ENABLE . HANDSH 
AKE) 

66  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

67  CALL 

I . V ( I , DISABLE . HANDSHAKE , SET . THRESHOLD , SET . START . BIT , SET . NUM . 
BITS) 

68  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

69  CALL  I. V( I, SET. LOGIC. SENSE, J, J, J) 

70  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

71  CALL  I. C( I, POSITIVE, NEGATIVE, TWOS, UNSIGNED) 

72  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

73  CALL  I.C(I,OC,TTL,R0,Rl) 

74  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

75  CALL  I.C(I,R2,R3,R4,R5) 

76  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

77  CALL  I.C(I,R6,R7,R8,R9) 

78  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

79  CALL  I.C(I,R10,R11,R12,R13) 

80  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

81  CALL  I.C(I,R14,R15,R16,J) 

82  IF    PCIB.ERRO0    THEN   ERROR   PCIB.BASERR 

83  1=6 

84  CALL 

I . V ( I , SET . FREQUENCY , SET . AMPLITUDE , SET . OFFSET , SET . SYMMETRY ) 
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8  5  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

86  CALL  I. V(I, SET. BURST. COUNT, J, J, J) 

87  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

88  CALL  I. C( I, SINE, SQUARE, TRIANGLE, CONTINUOUS) 

89  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

90  CALL  I. C(I, GATED, BURST, J, J) 

91  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

92  1=7 

93  CALL 

I . V ( I , AUTOSCALE , CALIBRATE , SET . SENSITIVITY , SET . VERT . OFFSET) 

94  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

95  CALL 

I . V ( I , SET . COUPLING , SET . POLARITY , SET . SWEEPSPEED , SET . DELAY ) 

96  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

97  CALL 

I . V ( I , SET . TRIG . SOURCE , SET . TRIG . SLOPE , SET . TRIG . LEVEL , SET . TRIG 
.MODE) 

98  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

99  CALL 

I . V ( I , GET . SINGLE . WF , GET . TWO . WF , GET . VERT . INFO , GET . TIMEBASE . IN 
FO) 

100  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

101  CALL 

I . V ( I , GET . TRIG . INFO , CALC . WFVOLT , CALC . WFTIME , CALC . WF . STATS ) 

102  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

103  CALL 

I . V ( I , CALC . RISETIME , CALC . FALLTIME , CALC . PERIOD , CALC . FREQUENCY 

) 

104  IF    PCIB.ERROO    THEN   ERROR   PCIB.BASERR 

105  CALL 

I . V ( I , CALC . PLUSWIDTH , CALC . MINUS WIDTH , CALC . OVERSHOOT , CALC . PRE 
SHOOT) 

106  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

107  CALL 

I . V ( I , CALC . PK . TO . PK, SET . TIMEOUT , SCOPE . START , MEASURE . SINGLE . W 

F) 

108  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

109  CALL  I.V(I,MEASURE.TWO.WF,J,J,J) 

110  IF    PCIB.ERROO    THEN   ERROR   PCIB.BASERR 

111  CALL    I.C(I,R10NANO,R100NANO,R1MICRO,R10MICRO) 

112  IF    PCIB.ERROO    THEN   ERROR   PCIB.BASERR 

113  CALL    I.C(I,R100MICRO,R1MILLI,R10MILLI,R100MILLI) 

114  IF    PCIB.ERROO    THEN    ERROR   PCIB.BASERR 

115  CALL    I.C(I,R1,R10,R2  0NANO,R200NANO) 

116  IF   PCIB.ERROO   THEN   ERROR   PCIB.BASERR 

117  CALL    I.C(I,R2MICRO,R2  0MICRO,R200MICRO,R2MILLI) 

118  IF   PCIB.ERROO   THEN   ERROR   PCIB.BASERR 

119  CALL   I.C(I,R20MILLI,R200MILLI,R2,R20) 

120  IF    PCIB.ERROO    THEN   ERROR   PCIB.BASERR 

121  CALL    I.C(I,R50NANO,R500NANO,R5MICRO,R50MICRO) 

122  IF    PCIB.ERROO    THEN    ERROR   PCIB.BASERR 
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123  CALL  I.C(I,R500MICRO,R5MILLI,R50MILLI,R500MILLI) 

124  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

125  CALL  I.C(I,R5,R50,CHAN.A,CHAN.B) 

126  IF    PCIB.ERROO    THEN   ERROR   PCIB.BASERR 

127  CALL  I. C( I, EXTERNAL, POSITIVE, NEGATIVE, AC) 

128  IF    PCIB.ERROO   THEN   ERROR   PCIB.BASERR 

129  CALL  I. C( I, DC, TRIGGERED, AUTO. TRIG, AUTO. LEVEL) 

130  IF   PCIB.ERROO   THEN   ERROR   PCIB.BASERR 

131  CALL  I. C( I, XI, XI 0, STANDARD, AVERAGE) 

132  IF   PCIB.ERROO   THEN   ERROR   PCIB.BASERR 

133  1=8 

134  CALL  I. V( I, OPEN. CHANNEL, CLOSE. CHANNEL, J, J) 

135  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
13  6  CALL  C.S 

137  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

138  I$=PCIB.DIR$+"\PCIB.PLD" 

139  CALL  L.P(I$) 

140  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

141  I$="DMM.01":I=3:J=0:K=0:L=1 

142  CALL  DEFINE(DMM.01,I$,I,J,K,L) 

143  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

144  I$="Func.Gen.01":I=6sJ=0:K=l:L=l 

145  CALL  DEFINE(FUNC.GEN.01,I$,I,J,K,L) 

146  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

147  I$="Scope . 01" : 1=7 : J=0 : K=2 : L=l 

148  CALL  DEFINE(SCOPE.01,I$,I,J,K,L) 

149  IF   PCIB.ERROO   THEN   ERROR   PCIB.BASERR 

150  I$="Counter . 01" : 1=1 : J=0 : K=3 : L=l 

151  CALL  DEFINE (COUNTER. 01,1$, I, J, K,L) 

152  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

153  I$="Dig. In. 01" : 1=4 : J=0 : K=4 : L=l 

154  CALL  DEFINE(DIG.IN.01,I$,I,J,K,L) 

155  IF   PCIB.ERROO   THEN   ERROR   PCIB.BASERR 

156  I$="Dig. Out . 01" : 1=4 : J=l : K=4 : L=l 

157  CALL  DEFINE(DIG.OUT.01,I$,I,J,K,L) 

158  IF   PCIB.ERROO   THEN   ERROR   PCIB.BASERR 

159  I$="Relay . Act . 01" : 1=8 : J=0 : K=5 : L=l 

160  CALL  DEFINE (RELAY. ACT. 01,1$, I, J, K,L) 

161  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

162  I$="Relay.Mux.01":I=2:J=0:K=6:L=l 

163  CALL  DEFINE (RELAY. MUX. 01,1$, I, J, K,L) 

164  IF   PCIB.ERROO   THEN   ERROR   PCIB.BASERR 

800  I$=ENVIR0N$( "PANELS ")+"\PANELS. EXE" 

801  CALL  L.S(I$) 

899  GOTO  2 

900  IF  ERR=PCIB.BASERR  THEN  GOTO  903 

901  PRINT  "BASIC  error  #";ERR;"  occurred  in  line  ";ERL 

902  STOP 

903  TMPERR=PCIB.ERR:IF  TMPERR=0  THEN  TMPERR=PCIB. GLBERR 

904  PRINT  "PC  Instrument  error  #";TMPERR;"  detected  at  line 
"  ;  ERL 
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905  PRINT  "Error:  ";PCIB.ERR$ 

9  06  IF  LEFT$(PCIB.NAME$/1)<>CHR$(32)  THEN  PRINT  "Instrument: 

";PCIB.NAME$ 

907  STOP 

908  COMMON  PCIB. DIR$ , PCIB. SEG 

909  COMMON  LD. FILE, GET. MEM, PANELS , DEF. ERR 

910  COMMON 

PCIB . BASERR , PCIB . ERR, PCIB . ERR$ , PCIB . NAME$ , PCIB . GLBERR 

911  COMMON 

READ . REGISTER , READ . SELFID , DEFINE , INITIALIZE . SYSTEM , ENABLE . SY 
STEM , DISABLE . SYSTEM , INITIALIZE , POWER . ON , MEASURE , OUTPUT , START 
, HALT , ENABLE . INT . TRIGGER , DISABLE . INT . TRIGGER , ENABLE . OUTPUT , D 
I SABLE . OUTPUT , CHECK . DONE , GET . STATUS 

912  COMMON 

SET . FUNCTION , SET . RANGE , SET . MODE , WRITE . CAL , READ . CAL , STORE . CAL 
, DELAY , SAVE . SYSTEM , SET . GATETIME , SET . SAMPLES , SET . SLOPE , SET . SO 
URCE , ZERO . OHMS , SET . SPEED , SET . COMPLEMENT , SET . DRIVER , OUTPUT . NO 
.WAIT , ENABLE . HANDSHAKE , DISABLE . HANDSHAKE 

913  COMMON 

SET . THRESHOLD , SET . START . BIT , SET . NUM . BITS , SET . LOGIC . SENSE , SET 
. FREQUENCY , SET . AMPLITUDE , SET . OFFSET , SET . SYMMETRY , SET . BURST . C 
OUNT , AUTOSCALE , CALIBRATE , SET . SENSITIVITY , SET . VERT . OFFSET , SET 
. COUPLING , SET . POLARITY , SET . SWEEPSPEED 

914  COMMON 

SET . DELAY , SET . TRIG . SOURCE , SET . TRIG . SLOPE , SET . TRIG . LEVEL , SET . 
TRIG . MODE , GET . SINGLE . WF , GET . TWO . WF , GET . VERT . INFO , GET . TIMEBAS 
E . INFO , GET . TRIG . INFO , CALC . WFVOLT , CALC . WFTIME , CALC . WF . STATS , C 
ALC . RISETIME , CALC . FALLTIME , CALC . PERIOD 

915  COMMON 

CALC . FREQUENCY , CALC . PLUSWIDTH , CALC . MINUSWIDTH , CALC . OVERSHOOT 
, CALC . PRESHOOT , CALC . PK . TO . PK, SET . TIMEOUT , SCOPE . START , MEASURE 
. SINGLE . WF , MEASURE . TWO . WF , OPEN . CHANNEL , CLOSE . CHANNEL 

916  COMMON 

FREQUENCY , AUTO . FREQ , PERIOD , AUTO . PER, INTERVAL, RATIO , TOTALIZE , 
R100MILLI , Rl , RIO , R100 , R1KILO , Rl OMEGA, RIO OMEGA, CHAN . A, CHAN . B , 
POSITIVE , NEGATIVE , COMN , SEPARATE , DCVOLTS , ACVOLTS , OHMS , R2  0  OMIL 
LI , R2 , R2  0 , R2  00 , R2KILO , R2  OKILO , R2  00KILO 

917  COMMON 

R2MEGA , R2  OMEGA , AUTOM , R2 . 5 , R12 . 5 , POSITIVE , NEGATIVE , TWOS , UNSIG 
NED , OC , TTL, RO , Rl , R2 , R3 , R4 , R5 , R6 , R7 , R8 , R9 , RIO , Rll , R12 , R13 , R14 
, R15 , R16 , SINE , SQUARE , TRIANGLE , CONTINUOUS , GATED , BURST , Rl ONANO 
, R100NANO , R1MICRO , R10MICRO , R100MICRO 

918  COMMON 

R1MILLI , R10MILLI , R100MILLI , Rl , RIO , R2  ONANO , R2  00NANO , R2MICRO , R 
2  OMICRO , R2  0  OMICRO , R2MI LLI , R2  OMI LLI , R2  0  OMI LLI , R2 , R2  0 , R5  ONANO , 
R500NANO , R5MICRO , R50MICRO , R500MICRO , R5MILLI , R50MILLI , R500MIL 
LI , R5 , R50 , CHAN . A , CHAN . B , EXTERNAL, POSITIVE 

919  COMMON 

NEGATIVE , AC , DC , TRIGGERED , AUTO . TRIG , AUTO . LEVEL, XI , XI 0 , STANDAR 
D , AVERAGE 
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SET  HP  BOXES  TO  INITIAL  SETTINGS  FOR  DATA  ACQUISITION 


92  0  COMMON 

DMM. 01 ,FUNC. GEN. 01, SCOPE. 01, COUNTER. 01, DIG. IN. 01, DIG. OUT. 01, 

RELAY . ACT .01, RELAY . MUX . 0 1 

999  'End  PCIB  Program  Shell 

1000  '  PROGRAM  "CALIB"RATE 
1010  OPTION  BASE  1 

1020  WANTAOA  =  0 

1030  ' 

1040  ' 

1050  ' 

1060  ' 

1070  ' 

1080  FILE$  =  "HOTWIRE. HPC" 

1090  CALL  INITIALIZE. SYSTEM (FILE$) 

1100  CALL  ENABLE . SYSTEM 

1110  CALL  SET. FUNCTION (DMM. 01, DCVOLTS) 

1120  CALL  SET.RANGE(DMM.'01,AUTOM) 

1130  CALL  SET.SPEED(DMM.01,R2.5) 

1140  CALL  ENABLE. INT. TRIGGER (DMM. 01) 

1141  CALL  SET.NUM.BITS(DIG.IN.01,R16) 

1142  LEVEL  =  5 

1143  CALL  SET.THRESHOLD(DIG.IN.01,LEVEL) 
1150 
1160 
1170 
1180 
1190 

1200  BEEP 
1210  CLS 

1220  COLOR  4:  LOCATE  5,14:  PRINT 
******************************************************** 

1230  LOCATE  9,14:  PRINT  "******************** 
********************** 

1240  COLOR  1:  LOCATE  9,37:  PRINT  "PART  1" 

COLOR  4:  PRINT  "**************" 

1250  LOCATE  6,14:  PRINT  "*":  LOCATE  7,14: 


PRINT  FIRST  SCREEN 


LOCATE  10,33: 


PRINT  "*":  LOCATE 


PRINT  "*":  LOCATE  7,65:  PRINT  "*":  LOCATE 


8,14:  PRINT  "*" 

12  60  LOCATE  6,65 

8,65:  PRINT  "*" 

1270  COLOR  1:  LOCATE  7,18:  PRINT  "WELCOME  TO  THE  HOT  WIRE 

FLOW  MEASUREMENT  LAB" 

1280  LOCATE  12,27:  PRINT  "THESIS  BY  LT  GREG  DOREMUS" 

1290  LOCATE  15,20:  PRINT  "THIS  PROGRAM  IS  USED  IN 

CALIBRATING  THE  " 

1300  LOCATE  17,13:  PRINT  "HOT  WIRE  SYSTEM  FOR  WIND  TUNNEL 

AIR  FLOW  MEASUREMENTS" 

1310  COLOR  4:  LOCATE  21,24 

1320  INPUT  "ENTER  1  TO  CONTINUE,  2  TO  ABORT" ;ANS1 

1330  IF  ANSI  <>  2  GOTO  1340  ELSE  5950 

1340  CLS:  LOCATE  10,36:  COLOR  20:  PRINT  "WAIT...":  COLOR  1 
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PRINT  SECOND  SCREEN 


1350  LOCATE  14,22:  PRINT  "AIRFOIL  IS  BEING  POSITIONED  TO 

ZERO" 

1360  LOCATE  16,20:  PRINT  "DEGREES  AOA  FOR  CALIBRATION 

MEASUREMENTS" 

1370  GOSUB  4940 

1380  GOSUB  5070 

1390 

1400 

1410 

1420 

1430 

1440  CLS:  LOCATE  6,10:  PRINT  "TO  START  THIS  EXPERIMENT  YOU 

MUST  FIRST  SET  UP  THE  HOT  WIRE" 

1450  LOCATE  8,12:  PRINT  "SYSTEM.  THE  HOT  WIRE  LABORATORY 

MANUAL  WILL  EXPLAIN  THIS" 

1460  LOCATE  10,6:  PRINT  "PROCEDURE  IN  DETAIL.  ONCE  THE 

SYSTEM  SET  UP  IS  TO  YOUR  SATISFACTION" 

1470  LOCATE  12,10:  PRINT  "YOU  MUST  MANUALLY  START  THE  WIND 

TUNNEL  AND  SET  THE  SLOWEST" 

1480  LOCATE  14,16:  PRINT  "TUNNEL  VELOCITY  YOU  WILL  USE  IN 

YOUR  EXPERIMENT." 

1490  COLOR  4:  LOCATE  20,24:  PRINT  "ENTER  1  TO  CONTINUE,  2  TO 

ABORT" 

1500  LOCATE  21,28:  INPUT  "OR  3  TO  GO  BACK  ONE  PAGE";ANS2 

1510  COLOR  1:  ON  ANS2  GOTO  1560,5950,1210 

1520 

1530 

1540 

1550 

1560  CLS:  LOCATE  10,22:  PRINT  "IF  HOTWIRE  SYSTEM  IS  SET  UP 

PROPERLY" 

1570  COLOR  4:  LOCATE  12,29:  PRINT  "START  WIND  TUNNEL  NOW" 

1580  LOCATE  20,24:  PRINT  "ENTER  1  TO  CONTINUE,  2  TO  ABORT" 

1590  LOCATE  21,28:  INPUT  "OR  3  TO  GO  BACK  ONE  PAGE";  ANS3 

1600  COLOR  1:  ON  ANS3  GOTO  1610,5950,1440 

1610 

1620 

1630 

1640 

1650 

1660  CLS:  LOCATE  8,13:  PRINT  "A  SERIES  OF  WIND  VELOCITY 

MEASUREMENTS  IS  TO  BE  TAKEN" 

1670  LOCATE  10,12:  PRINT  "AT  WIND  SPEEDS  VARYING  FROM  THE 

SLOWEST  TUNNEL  VELOCITY" 

1680  LOCATE  12,11:  PRINT  "TO  BE  USED  IN  YOUR  EXPERIMENT  UP 

TO  THE  HIGHEST  VELOCITY." 

1690  LOCATE  14,15:  PRINT  "A  MINIMUM  OF  10  READINGS  SHOULD  BE 

TAKEN  TO  CREATE" 

1700  LOCATE  16,29:  PRINT  "THE  CALIBRATION  CURVE." 

1710  COLOR  4:  LOCATE  2  0,24:  PRINT  "ENTER  1  TO  CONTINUE,  2  TO 

ABORT" 


PRINT  THIRD  SCREEN  -  START  WIND  TUNNEL 


PRINT  FOURTH  SCREEN 
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OPEN  DATA  FILES 


1720  LOCATE  21,28:  INPUT  "OR  3  TO  GO  BACK  ONE  PAGE";ANS4 

1730  COLOR  1:  ON  ANS4  GOTO  1790,5950,1560 

1740 

1750 

1760  '  ENTER  DATA  POINTS  FOR  CALIBRATION  CURVE 

1770 

1780 

1790  CLS:  LOCATE  12,12:  PRINT  "ENTER  THE 

YOU  WILL  TAKE  TO  CREATE" 

1800  COLOR  4:  LOCATE  12,22:  PRINT  "NUMBER  OF  DATA  POINTS" 

COLOR  1 

1810  LOCATE  14,19:  INPUT  "THE  CALIBRATION  CURVE  FOR  THIS 

EXPERIMENT. " ; POINTS 

182  0  CLS 

1830  TEST  =  1 

1840  CALL  OUTPUT (RELAY. MUX. 01, TEST) 

1850  DIM  Z(POINTS,2),  X ( POINTS ) ,  Y(POINTS),  C(10),  D(10), 

KNOTS (POINTS, 2) ,  DEC(16) 

1860 

1870 

1880 

1890 

1900 

1910  TIME$  =  "00:00:00" 

1920  COLOR  4 

1930  LOCATE  14,18:  PRINT  "CREATING  DATA  FILES  FOR 

CALIBRATION  VOLTAGES" 

1940  COLOR  1 

1950  OPEN  "HOTWIRE.DAT"  FOR  OUTPUT  AS  #1 

1960  CLOSE  #1 

1970  TIME  =  VAL(RIGHT$(TIME$,2) ) 

1980  IF  TIME  <  5  GOTO  1970  ELSE  1990 

1990 

2000 

2010 

2020 

2030 

2040  READING  =  0:  CHANNEL  =  0 

2050  FOR  READING  =  1  TO  POINTS 

2051  CLS:  LOCATE  10,15:  PRINT  "DEPRESS  BUTTON     ON  THE 
PARASCIENTIFIC  TRANSDUCER" 

2052  COLOR  4:  LOCATE  10,30:  PRINT  "PI":  COLOR  1 

2053  LOCATE  12,15:  PRINT  "MAKE  SURE  NO  OTHER  FUNCTION 
BUTTONS  ARE  DEPRESSED" 

2054  '  COLOR  4:  LOCATE  21,28:  PRINT  "PRESS  ENTER  TO 
CONTINUE" 

2055  A$  =  INKEY$:  IF  A$  =  ""  THEN  2055 
2060    CLS:  LOCATE  8,17:  COLOR  1 

2070    IF  READING  =  1  GOTO  2080  ELSE  2100 

2080    PRINT  "IS  THE  PROBE  POSITIONED  FOR  YOUR  FIRST 

READING?" 
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START  DATA  ACQUISITION 


GOTO  2110 


2090 

2100    LOCATE 

READING?" 


2110 

2120 

2130 

2140 

FOR  A 

2150 

2160 

2170 

2180 

WIRE 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320  • 

2330 

2340 

2342 

2344 

2370 

2410 

2420 

2430 

2440 

2442 

2450 

2460 

2470 

2480 

2490 

2500 

DATA" 

2510 

2520 

COLOR 

2530 

COLOR 

2540 


8,19:  PRINT  "IS  THE  PROBE  POSITIONED  FOR  A  NEW 
LOCATE  9,36:  INPUT  "1  =  YES";ANS5 


COLOR  4 

IF  ANS5  <>  1  GOTO  2060 
COLOR  1 

CLS:  LOCATE  16,19: 
NEW  FLOW  VELOCITY?" 
COLOR  4:  LOCATE 

IF  ANS6  <>  1  GOTO  2140 
CLS:  COLOR  20:  LOCATE  8 


COLOR  1:  PRINT  "IS  THE  TUNNEL  SET 


17,36:  INPUT 


"1  =  YES";ANS6 


35:  PRINT  "STAND  BY" 


COLOR  1:  LOCATE  13,24:  PRINT  "DATA  VALUES  FOR  HOT 
VOLTAGE" 

LOCATE  15,25 


LOCATE  17,30: 


PRINT 
PRINT 


"AND 
"NOW 


TEST  SECTION  VELOCITY  ARE" 
BEING  COLLECTED" 


FOR  CHANNEL  =  1  TO  2 

ON  CHANNEL  GOTO  2230,2342 

CALL  OUTPUT (RELAY. MUX. 01, CHANNEL) 
CALL  SET.RANGE(DMM.01,R20) 

1=0:  SUM1  =  0:  FOR  I  =  1  TO  10 
IF  I  =  1  GOTO  2270  ELSE  2290 
TIME  =  .5 
CALL  DELAY (TIME) 
CALL  MEASURE (DMM. 01, C (I) ) 
SUM1  =  SUM1  +  ABS(C(I)) 
NEXT  I 
AVG1  =  SUM1/10 
Z (READING, CHANNEL)  =  AVG1 
GOTO  2780 
1  =  0 
SUM2  =  0 
FOR  I  =  1  TO  10 

CALL  MEASURE(DIG.IN.01,D(I) ) 
SUM2  =  SUM2  +  D(I) 
NEXT  I 
DN  =  SUM2/10 
GOSUB  5480 

Z (READING, CHANNEL)  =  PRESS 
GOSUB  4270 
GOSUB  4700 

Z (READING, CHANNEL)  =  VEL 
KNOTS (READING, CHANNEL)  =  KTS 
BEEP:  CLS:  LOCATE  2,32:  PRINT  "CALIBRATION 


LOCATE 
LOCATE 


4,17: 
4,46: 


PRINT  "HOT  WIRE  VOLTAGE" 
PRINT  "WIND  TUNNEL  VELOCITY" 


LOCATE  6,19:  PRINT  USING  "#.####" ;Z (READING, 1) 
LOCATE  6,26:  PRINT  "VOLTS":  COLOR  4 
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2550  LOCATE  6,41:  PRINT  USING 

"###. ###" ;Z (READING, 2) :  COLOR  1 

2560  LOCATE  6,49:  PRINT  "FT/SEC":  COLOR  4 

2570  LOCATE  6,58:  PRINT  USING 

"###.###"; KNOTS (READING , 2 ) :  COLOR  1 

2580  LOCATE  6,66:  PRINT  "KNOTS":  LOCATE  6,56:  PRINT 

ii —ii 

2590 

"; READING;" 

2600 

2610 

2620 

COLOR  1 

2630 

2640 

2650 


LOCATE  8,15:  PRINT  "THIS  IS  DATA  POINT 
OF  "? POINTS 
LOCATE  10,32: 


LOCATE  12,15: 


PRINT  "WIND  TUNNEL  DATA" 
PRINT  "TUNNEL  TEMP  =" 


COLOR  4:  LOCATE  12,29:  PRINT  USING  "##.##";TF: 


LOCATE  12,35:  PRINT  "DEG  F" 

LOCATE  14,15:  PRINT  "TUNNEL  STATIC  PRESSURE  =' 
COLOR  4:  LOCATE  14,40:  PRINT  USING 
"####. ##";PSTAT:  COLOR  1 

2660  LOCATE  14,48:  PRINT  "LB/FTA2" 

2670  LOCATE  16,15:  PRINT  "TUNNEL  AIR  DENSITY  =" 

2680  COLOR  4:  LOCATE  16,36:  PRINT  USING 

"#.#####" ; DENSITY:  COLOR  1 


2690 

2700 

2710 

COLOR  1 

2720 

2730 

2740 

VELOCITY 

2750 

CONTINUE" 

2760 


IN 


LOCATE  16,44:  PRINT  "LB/FTA3" 

LOCATE  18,15:  PRINT  "DYNAMIC  PRESSURE  =" 

COLOR  4:  LOCATE  18,34:  PRINT  USING  "##.##" ;Q: 

LOCATE  18,40:  PRINT  "LB/FTA2" 
IF  READING  =  POINTS  GOTO  2750 
COLOR  20:  LOCATE  21,24:  PRINT  "SET  NEW  FLOW 
TUNNEL" 
COLOR  4:  LOCATE  23,27:  PRINT  "PRESS  ANY  KEY  TO 


A$  =  INKEY$:  IF  A$  = 
COLOR  1 
NEXT  CHANNEL 
NEXT  READING 
BNUM1  =  Z(l,2) 

Z (POINTS, 2) 
KNOTS (1,2) 
KNOTS < POINTS, 2) 


ii  ii 


THEN  2760 


BNUM2  = 
BNUM3  = 
BNUM4  = 


2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890  CLS:  LOCATE  8,15:  PRINT 

IN  DATA  FILE  LABELED" 

2900  COLOR  4:  LOCATE  10,34 

2910  I  =  0:  CI  =  0:  C2  =  0 

2920  FOR  I  =  1  TO  POINTS 

2930    X(I)  =  SQR(Z(I,2) ) 

2940    Y(I)  =  (Z(I,1))A2 


STORE  DATA  IN  DATA  FILE  LABELED  "HOTWIRE.DAT" 


"DATA  VALUES  ARE  BEING  STORED 


PRINT  "HOTWIRE.DAT": 
Dl  =  0:  D2  =  0 


COLOR  1 
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2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

YOUR 

3200 

3210 

3220 

3230 

3240 

3250 

COLOR 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

COPY 

3370 


CI  =  CI  +  X(I) 

C2  =  C2  +  X(I)A2 

Dl  =  Dl  +  Y(I) 

D2  =  D2  +  X(I)*Y(I) 
NEXT  I 

OPEN  "HOTWIRE.DAT"  FOR  APPEND  AS 
J  =  0:  FOR  J  =  1  TO  POINTS 

PRINT  #1,USING»###.####»;Y(J) , 
NEXT  J 
CLOSE  #1 


#1 
X(J) 


SOLVE  FOR  SLOPE  AND  Y-INTERCEPT  OF  REGRESSION  LINE 


A  =  (P0INTS*D2-C1*D1)/(P0INTS*C2-C1A2) 

IF  ABS(A)  <  1E-08  THEN  A=0 

B  =  (D1-A*C1) /POINTS 

IF  ABS(B)  <  1E-08  THEN  B=0 


PRINT  SUMMARY  OF  CALIBRATION  DATA  ON  SCREEN 


LOCATE  16,12:  PRINT  "WOULD  YOU  LIKE  TO  SEE  A  SUMMARY  OF 
CALIBRATION  DATA" 

LOCATE  18,33:  PRINT  "ON  THE  SCREEN?":  COLOR  4 
LOCATE  20,27:  INPUT  "ENTER  1  FOR  YES,  0  FOR  NO";ANS7 
IF  ANS7  <>  0  GOTO  3230  ELSE  3640 

COLOR  1:  CLS:  LOCATE  2,32:  PRINT  "CALIBRATION  DATA" 
LOCATE  4,9:  PRINT  "HOT  WIRE  VOLTAGE  (VOLTS)" 
LOCATE  4,41:  PRINT  "WIND  TUNNEL  VELOCITY  (FT/SEC)": 
4 


J  =  0 

FOR  J  =  1  TO  POINTS 
LOCATE  J+5,18:  PRINT 
LOCATE  J+5,52:  PRINT 

NEXT  J 


USING"#.####";Z(J,1) 
USING"###.###»;Z(J,2) 


SEND  CALIBRATION  DATA  TO  HP  LASER  JET  PRINTER 


PRINT  "WOULD  YOU  LIKE  A  HARD 


INPUT  "ENTER  1  FOR  YES,  2  FOR 


COLOR  1:  LOCATE  22,20: 

OF  THIS  DATA?" 

COLOR  4:  LOCATE  23,27: 
NO";ANS8 

3380  IF  ANS8  <>  2  GOTO  3390  ELSE  3640 

3  390  CLS:  LOCATE  12,26:  PRINT  "DATA  SENT  TO  LASER  PRINTER" 
COLOR  1 
3400  TIME$  =  "00:00:00" 
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TIME  =  VAL(RIGHT$(TIME$/2) ) 

IF  TIME  <  5  GOTO  3410  ELSE  3430 

LPRINT:  LPRINT 

LPRINT  TAB (32):  LPRINT  "CALIBRATION  DATA" 

LPRINT 

LPRINT  TAB (9):  LPRINT  "HOTWIRE  VOLTAGE 

TUNNEL  VELOCITY" 

LPRINT  TAB(42):  LPRINT  "(FT/SEC) 

LPRINT 

J  =  0 

FOR  J  =  1  TO  POINTS 

LPRINT  TAB (18) ; 

LPRINT  USING"#.####";Z(J,1) ; 

LPRINT  TAB (43) ; 

LPRINT  USING"###.###";Z(J,2) ; 

LPRINT  TAB (61) ; 

LPRINT  USING" ###:###", -KNOTS (J, 2) 
NEXT  J 
LPRINT  CHR$(12) 


3410 

3420 

3430 

3440 

3450 

3460 

WIND 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3  640  COLOR 

SEE  A  PLOT 

3  650  COLOR 

NO";ANS9 

3  660  COLOR 


(VOLTS) 

(KTS) " 


PRINT  VIEW  PLOT  OF  CALIBRATION  CURVE  ON  THE  SCREEN 


1:  CLS:  LOCATE  12,13:  PRINT  "WOULD  YOU  LIKE  TO 

OF  THE  CALIBRATION  CURVE?" 

4:  LOCATE  14,27:  INPUT  "ENTER  1  FOR  YES,  2  FOR 


CLS: 
i 


1:  IF  ANS9  <>  2  GOTO  3  670  ELSE  3730 
SHELL  "COMMAND/C  HOT1.BAT" 


3670 

3680 

3690 

3700 

PRINTER 

3710  ' 

3720  ■ 

373  0  COLOR  1:  CLS:  LOCATE 
HARD  COPY  OF  THIS  GRAPH?" 

374  0  LOCATE  14,27:  COLOR  4 
NO";ANS10 

3750  COLOR  1:  IF  ANS10  <> 
3760  CLS:  COLOR  1:  LOCATE 
AUTOMATICALLY  LEAVE" 

3770  LOCATE  12,22:  PRINT  "BASICA  AND  ENTER  A  GRAPHICS 

PROGRAM." 

3780  LOCATE  14,14:  PRINT  "THE  PLOT  TAKES  APPROXIMATELY  4-6 

MINUTES  TO  COMPLETE" 

3790  COLOR  4:  LOCATE  21,27:  PRINT  "PRESS  ANY  KEY  TO 

CONTINUE" 

3800  A$  =  INKEY$:  IF  A$  =  ""  THEN  3800 

3810  COLOR  1 


SEND  HARD  COPY  OF  CALIBRATION  CURVE  TO  HP  LASER  JET 


12,19:  PRINT  -"WOULD  YOU  LIKE  A 

:  INPUT  "ENTER  1  FOR  YES,  2  FOR 

2  GOTO  3760  ELSE  3840 

10,21:  PRINT  "THIS  PROGRAM  WILL 
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3820  SHELL  "HOT2.BAT" 

3830  BEEP 

384  0  OPEN  "CHAIN. DAT"  FOR  OUTPUT  AS  #1 

3850    PRINT  #1,USING"####.###";  A,  B, 

BNUM4,  TF 

3860  CLOSE  #1 

3870 

3880 

3890 

RANGE 

3900 

3910 

3920  CLS:  LOCATE  8,14: 

CURVE  PORTION  OF  YOUR" 

393  0  LOCATE  10,10:  PRINT  "EXPERIMENT. 

NOW  ACCURATELY  CALIBRATED" 

3940  LOCATE  12,13:  PRINT  "WITHIN  A  WIND  TUNNEL  VELOCITY 

RANGE  OF" 

3950  COLOR  4:  LOCATE  12,52:  PRINT  USING"### . ###" ;Z (1, 2) 
12,61:  PRINT  "FT/SEC":  COLOR  1 
14 , 17 :  PRINT  " ( " :  LOCATE  14,18:  PRINT  USING 


BNUM1 ,  BNUM2 ,  BNUM3 , 


PRINT  CONCLUDING  SCREEN  AND  SUMMARY  OF  CALIBRATION 


PRINT  "THIS  CONCLUDES  THE  CALIBRATION 


THE  HOT  WIRE  SYSTEM  IS 


3960  LOCATE 
3970  LOCATE 


"###.###";KNOTS(l,2) 
3980  LOCATE  14,26:  PRINT  "KTS) 
3990  LOCATE  14,34:  PRINT  USING 
14,42:  PRINT  "FT/SEC":  COLOR  1: 
4000  LOCATE  14,50:  PRINT  USING 
4010  LOCATE  14,58:  PRINT  "KTS)" 
402  0  COLOR  4:  LOCATE  20,24:  PRINT 
ABORT" 

4  03  0  LOCATE  21,28:  INPUT  "OR  3  TO 
ON  ANS11  GOTO  4100,5950,3730 


TO":  COLOR  4 

»###. ###";Z (POINTS, 2) :  LOCATE 

LOCATE  14,49:  PRINT  "(" 
»•###.  ###»  ;KNOTS  (POINTS,  2) 


"ENTER  1  TO  CONTINUE,  2  TO 
GO  BACK  ONE  PAGE";ANS11 


PRINT  LAST  SCREEN  AND  CALIBRATION  EQUATION 


4040 

4050  ' 

4060  ' 

4070  ' 

4080  ■ 

4090  ' 

4100  CLS:  LOCATE  8,30:  PRINT  "Y  =  " 

4110  LOCATE  8,34:  PRINT  USING  "##.###" ;A:  LOCATE 

"X" 

4120  IF  B>=0  THEN  LOCATE  8,42:  PRINT  "+" 

USING  »##.###»;B  ELSE  LOCATE  8,42:  PRINT 

PRINT  USING  "###.###" ;ABS(B) 

413  0  COLOR  1:  LOCATE  10,10:  PRINT  "THIS  IS 

EQUATION  OF  THE  CALIBRATION  CURVE" 

4140  LOCATE  12,19:  PRINT  "IT  WILL  AUTOMATICALLY  BE  ENTERED 

IN  PART  2" 

4150  LOCATE  14,26:  PRINT  "OF  YOUR 

4160  COLOR  4:  LOCATE  2  0,27:  PRINT 

TWO" 

4170  LOCATE  21,32:  PRINT  "ENTER  2  TO  ABORT" 


8,40 


LOCATE  8,44: 
"-":  LOCATE  8 


PRINT 

PRINT 
44: 


THE  STRAIGHT  LINE 


HOT  WIRE  EXPERIMENT" 

"  ENTER  1  TO  START  PART 


75 


4180 

PAGE" 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 


LOCATE  22,25:  INPUT  "OR  ENTER  3  TO  GO  BACK  ONE 

;ANS12 

COLOR  1 

ON  ANS12  GOTO  4210,5950,3920 

LOAD  "PART2",R 


SUBROUTINE  FOR  MEASURING  WIND  TUNNEL  TEMPERATURE 


Al  =  38.709457# 

A2  =  .037085566# 

A3  =  5.649552E-05 

BO  =  .10086091# 

Bl  =  25727. 94369# 

B2  =  -767345. 8295# 

B3  =  78025595. 81# 

B4  =  -9247486589# 

B5  =  6.97688E+11 

B6  =  -2.66192E+13 

B7  =  3.94078E+14 

B8  =  0 

B9  =  0 

REF.OUT  =  8 

CALL  OUTPUT (RELAY. MUX. 01, REF.OUT) 

CALL  ENABLE . OUTPUT (RELAY . MUX .01) 

CALL  SET. FUNCTION (DMM. 01, DCVOLTS) 

CALL  SET.RANGE(DMM.01,R2) 

CALL  DISABLE .JNT. TRIGGER (DMM. 01) 

COUNT  =  .5 

TIME$  -  "00:00:00" 

IF  TIMER  <  COUNT  GOTO  4480  ELSE  4490 
CALL  MEASURE (DMM. 01, V) 
TR  =  100*V 

ER  =  (TR*(A1+TR*(A2+TR*A3) ) )*10A-6 
CALL  SET. RANGE (DMM. 01,R200MILLI) 
TC  =  7 
CALL  OUTPUT (RELAY. MUX. 01, TC) 

COUNT  =  .5 

TIME$  =  "00:00:00" 

IF  TIMER  <  COUNT  GOTO  4570  ELSE  4580 
CALL  MEASURE (DMM. 01, ET) 
E  =  ER  +  ET 

Z  =  B5+E*(B6+E*(B7+E*(B8+E*B9) ) ) 
T  =  B0+E*(B1+E*(B2+E*(B3+E*(B4+E*Z) ) ) ) 
TF  =  (1.8*T)  +■  32 
RETURN 


•  SUBROUTINE  FOR  WIND  TUNNEL  DATA: 

'  STATIC  PRESS  IN  TEST  SECTION 


76 


AIR  DENSITY  IN  TEST  SECTION 
DYNAMIC  PRESSURE  IN  TEST  SECTION 
TEST  SECTION  VELOCITY,  FT/SEC  AND 


4680 

4690 

4700 

4710 

4720 

4722  CLS:  COLOR  1:  LOCATE  10,15 

THE  PARASCIENTIFIC  TRANSDUCER" 

4724  COLOR  4:  LOCATE  10,30:  PRINT  "P2" 

4726  LOCATE  12,15:  PRINT  "MAKE  SURE  NO 

BUTTONS  ARE  DEPRESSED" 

4728  COLOR  4:  LOCATE  21,28:  PRINT  "PRESS 

A$  =  INKEY$:  IF  A$  =  ""  THEN  4730 

CALL  MEASURE (DIG. IN. 01, DN) 

GOSUB  5480 

PSTAT  =  PRESS* 14 4 

TTEMP  =  TF  +  459.67 

GASCONST  =  53.3 

DENSITY  =  PSTAT/ (GASCONST*TTEMP) 

DELTAP  =  Z (READING, CHANNEL) * 14 4 

Q  =  DELTAP/. 9 3 

IF  Q  <  1  THEN  Q  =  0 

GC  =  32.174 

VEL  =  SQR(ABS( (2*Q*GC) /DENSITY) ) 

KTS  =  VEL  *  .5924 

RETURN 


KTS 


PRINT  "DEPRESS  BUTTON 


COLOR  1 
OTHER  FUNCTION 


ON 


ENTER  TO  CONTINUE" 


4730 
4800 
4802 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
5110 
5120 
5130 
5140 
5150 
5160 
5170 


SUBROUTINE  FOR  READING  AOA 


INPUT. CH  =  5 

CALL  OUTPUT (RELAY. MUX. 01, INPUT. CH) 

CALL  ENABLE. OUTPUT (RELAY. MUX. 01) 

TIME  =  3 

CALL  DELAY (TIME) 

CALL  MEASURE ( DMM .01, AOA ) 

HAVEAOA  =  (AOA-.003) *1000/10 

RETURN 


SUBROUTINE  FOR  AOA  SELECTION 


WHILE  (CINT(ABS(WANTAOA  -  HAVEAOA)  *  10)  /  10)  >=  .25 
ABSDIFF  =  ABS(WANTAOA  -  HAVEAOA) 
IF  ABSDIFF  >  36  THEN  ABSDIFF  =36 
SWITCH. 1  =  1 
SWITCH. 2  =  2 

CALL  OPEN. CHANNEL (RELAY . ACT . 01 , SWITCH . 1) 
CALL  OPEN . CHANNEL ( RELAY . ACT .01, SWITCH . 2 ) 
IF  WANTAOA  >  HAVEAOA  GOTO  5210  ELSE  5340 
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5180      ' 

5190      '  THIS  PART  OF  THE  SUBROUTINE  MOVES  THE  TURNTABLE 

FORWARD 

5200      ' 

5210      IF  (CINT(ABSDIFF*10)/10)  >=  .4  THEN  TIMING  =  100 

ELSE  TIMING  =  1 

5220      COUNT  =  ( (ABSDIFF  *  TIMING)  +  190)  /  1000 

5230      IF  COUNT  >  3  THEN  COUNT  =  3 

5240      CALL  CLOSE. CHANNEL (RELAY. ACT. 01, SWITCH. 1) 

5250      CALL  DELAY (COUNT) 

5260      CALL  OPEN. CHANNEL (RELAY. ACT. 01, SWITCH. 1) 

5270        TIME  =  .5 

5280        CALL  DELAY (TIME) 

5290      GOSUB  4940 

5300      GOTO  5440 

5310        ' 

5320        '  THIS  PART  OF  THE  SUBROUTINE  MOVES  THE  TURNTABLE 

IN  REVERSE 

5330        ■ 

5340        IF  (CINT(ABSDIFF*10)/10)  >=  .4  THEN  TIMING  =  100 

ELSE  TIMING  =  1 

5350        COUNT  =  ((ABSDIFF*  TIMING)  +  190)/1000 

5360        IF  COUNT  >  3  THEN  COUNT  =  3 

5370        CALL  CLOSE. CHANNEL (RELAY. ACT. 01, SWITCH. 2) 

5380        CALL  DELAY (COUNT) 

5390        CALL  OPEN. CHANNEL (RELAY. ACT. 01, SWITCH. 2) 

5400        TIME  =  .5 

5410        CALL  DELAY (TIME) 

5420        GOSUB  4940 

5430        GOTO  5440 

5440  WEND 

5450  RETURN 

5460  ■ 

5470  ■ 

5480  '  SUBROUTINE  FOR  CONVERTING  DECIMAL  PRESSURE  VALUES  TO 

NUMERIC  VALUES 

5490  .« 

5500  ' 

5510  COUNT  =  0:  PRESS  =  0 

5520  • 

553  0  '  DN  IS  THE  MEASURED  DECIMAL  NUMBER 

5540  '  DEC  IS  THE  DECIMAL  TRANSFORM  MATRIX 

5550  •  PRESS  IS  THE  CONVERTED  NUMERIC  PRESSURE  VALUE 

5560  ' 

5570  DEC(16)  =  10:  DEC (15)  =  8:  DEC (14)  =  4:  DEC (13)  =  2: 

DEC (12)  =  1 

5580  DEC(ll)  =  .8:  DEC(IO)  =  .4:  DEC(9)  =  .2:  DEC(8)  =  .1: 

DEC(7)  =  .08 

5590  DEC(6)  =  .04:  DEC(5)  =  .02:  DEC(4)  =  .01:  DEC(3)  = 

8.000001E-03 

5600  DEC(2)  =  .004:  DEC(l)  =  .002 
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5610  ' 

562  0  IF  DN  <  0  THEN  PRESS  =  DEC (16)  ELSE  GOTO  564  0 

5630  COUNT  =  -32768! 

5640  IF  (DN-COUNT)/16384  >=  1  THEN  PRESS  =  PRESS  +  DEC (15) 

ELSE  GOTO  5660 

5650  COUNT  =  COUNT  +  16384 

5660  IF  (DN-COUNT)/8192  >=  1  THEN  PRESS  =  PRESS  +  DEC (14) 

ELSE  GOTO  5680 

5670  COUNT  =  COUNT  +  8192 

5680  IF  (DN-COUNT)/4096  >=  1  THEN  PRESS  =  PRESS  +  DEC (13) 

ELSE  GOTO  5700 

5690  COUNT  =  COUNT  +  4096 

5700  IF  (DN-COUNT)/2048  >=  1  THEN  PRESS  =  PRESS  +  DEC (12) 

ELSE  GOTO  5720 

5710  COUNT  =  COUNT  +  2048 

5720  IF  (DN-COUNT)/1024  >=  1  THEN  PRESS  =  PRESS  +  DEC (11) 

ELSE  GOTO  5740 

5730  COUNT  =  COUNT  +  1024 

5740  IF  (DN-COUNT)/512  >=  1  THEN  PRESS  =  PRESS  +  DEC (10) 

ELSE  GOTO  5760 

5750  COUNT  =  COUNT  +  512 

5760  IF  (DN-COUNT)/256  >=  1  THEN  PRESS  =  PRESS  +  DEC (9)  ELSE 

GOTO  5780 

5770  COUNT  =  COUNT  +256 

5780  IF  (DN-COUNT)/128  >=  1  THEN  PRESS  =  PRESS  +  DEC (8)  ELSE 

GOTO  5800 

5790  COUNT  =  COUNT  +  128 

5800  IF  (DN-COUNTJ/64  >=  1  THEN  PRESS  =  PRESS  +  DEC (7)  ELSE 

GOTO  5820 

5810  COUNT  =  COUNT  +  64 

5820  IF  (DN-COUNT)/32  >=  1  THEN  PRESS  =  PRESS  +  DEC (6)  ELSE 

GOTO  5840 

583  0  COUNT  =  COUNT  +32 

5840  IF  (DN-COUNT)/16  >=  1  THEN  PRESS  =  PRESS  +  DEC (5)  ELSE 

GOTO  5860 

5850  COUNT  =  COUNT  +  16 

5860  IF  (DN-COUNT)/8  >=  1  THEN  PRESS  =  PRESS  +  DEC(4)  ELSE 

GOTO  5880 

587  0  COUNT  =  COUNT  +  8 

5880  IF  (DN-COUNT)/4  >=  1  THEN  PRESS  =  PRESS  +  DEC (3)  ELSE 

GOTO  5900 

5890  COUNT  =  COUNT  +  4 

5900  IF  (DN-COUNT)/2  >=  1  THEN  PRESS  =  PRESS  +  DEC (2)  ELSE 

GOTO  5920 

5910  COUNT  =  COUNT  +  2 

592  0  IF  (DN-COUNT)  =  1  THEN  PRESS  =  PRESS  +  DEC(l) 

593  0  RETURN 
5940  STOP 

5950  CLS:  SYSTEM 
5960  END 
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APPENDIX  C 


Program  Listing:  PART2 


1000  '  PROGRAM  "PART2" 

1010  CALL  DEF.ERR(PCIB.ERR,PCIB.ERR$,PCIB.NAME$,PCIB.GLBERR) 

1020  FILE$  =  " HOTWIRE. HPC" 

103  0  CALL  INITIALIZE. SYSTEM (FILE$) 

104  0  CALL  ENABLE . SYSTEM 

1050  CALL  SET. FUNCTION (DMM. 01 ,DCVOLTS) 

1060  CALL  SET.RANGE(DMM.01,AUTOM) 

1070  CALL  SET.SPEED(DMM.01,R2.5) 

1080  CALL  DISABLE. INT. TRIGGER (DMM. 01) 

109  0  CALL  ENABLE . OUTPUT ( RELAY. MUX. 01) 

1100  CALL  SET.NUM.BITS(DIG.IN.01,R16) 

1110  LEVEL  =  5 

1120  CALL  SET. THRESHOLD (DIG. IN. 01, LEVEL) 

1130  DIM  WANTAOA(50) ,  G(40),  HWVOLT(IOO),  POSITION(IOO) , 

V2(100),  X(100) ,Y(100) ,  PS.VEL(IOO),  S.PR(IOO),  Y2(100), 

SPRDATA(IOO) ,  DEC(16),  M(100),  Z(100) 

1140  DIM  YINT(IOO),  Y2INT(100),  F(70) 

CLS 

OPEN  "CHAIN. DAT"  FOR  INPUT  AS  #1 


1150 

1160 

1170 

1180 

1190 

1200 

F) 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 


INPUT  #1, 
CLOSE  #1 
TCAL  =  TF 
1  TS  IS  THE 

TS  =  482 


A,  B,  BNUM1,  BNUM2,  BNUM3 ,  BNUM4 ,  TF 


TEMPERATURE  OF  THE  HOTWIRE  SENSOR  (DEGREES 


PRINT  INTRO  SCREEN 


BEEP 
COLOR 


4:  LOCATE  5,14:  PRINT 


I'****************************************************" 
1290  LOCATE  9,14:  PRINT  "******************** 


a*******************" 

1300  COLOR  1:  LOCATE  9,37:  PRINT  "PART  2":  LOCATE  10,33: 


6,14 


***************** 

PRINT  "*":  LOCATE 


7,14:  PRINT  "*":  LOCATE 


COLOR  4:  PRINT 
1310  LOCATE 

:  PRINT  "*" 
LOCATE  6,65 
:  PRINT  "*" 
COLOR  1:  LOCATE  7,18 
MEASUREMENT  LAB" 
LOCATE  12,27:  PRINT  "THESIS  BY  LT  GREG  DOREMUS" 


8,14: 

1320 

8,65: 

1330 

FLOW 

1340 


PRINT  "*":  LOCATE  7,65:  PRINT  "*":  LOCATE 


PRINT  "WELCOME  TO  THE  HOT  WIRE 
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1350  LOCATE  15,20:  PRINT  "THIS  PROGRAM  IS  USED  TO  ANALYZE 

THE  WAKE" 

1360  LOCATE  16,18:  PRINT  "VELOCITY  PROFILES  OF  THE  NACA 

66(215)  -  216" 

1370  LOCATE  17,24:  PRINT  "AIRFOIL  USING  HOTWIRE  ANEMOMETRY" 

1380  COLOR  4:  LOCATE  21,24:  INPUT  "ENTER  1  TO  CONTINUE,  2  TO 

ABORT"; AN 1 

13  90  COLOR  1 

1400  IF  AN1  <>  2  GOTO  1440  ELSE  5180 

1410  ' 

1420  ■ 

143  0  '  INITIAL  AOA  SELECTION:  PREVIOUSLY  STORED  AOA  VALUES 

COME  FROM  DATA 

1440  '  FILE  LABELED  "AOA. DAT" 

1450  ' 

1460  ' 

1470  CLS:  LOCATE  9,9:  PRINT  "DO  YOU  WISH  TO  USE  THE  PREVIOUS 

ANGLE  OF  INCIDENCE  SELECTIONS?" 

1480  COLOR  4 

1490  LOCATE  20,27:  INPUT  "ENTER  1  FOR  YES,  0  FOR  NO" ;AN3 

1500  COLOR  1 

1510  IF  AN3  <>  0  GOTO  1520  ELSE  1740 

1520  LOCATE  13,28:  PRINT  "PREVIOUS  AOA  SELECTIONS:" 

1530  COLOR  4 

1540  OPEN  "AOA. DAT"  FOR  INPUT  AS  #1 

1550  SNUM  =  1 

1560  LOCATE  15,2 

1570  IF  EOF(l)  THEN  SNUM  =  -99 

1580  IF  SNUM  =  -99  GOTO  1630  ELSE  1590 

1590    INPUT  #1,  WANTAOA(SNUM) 

1600    IF  SNUM  =  1  THEN  PRINT  WANTAOA(SNUM) ;  ELSE  PRINT  ", 

";WANTAOA(SNUM) ; 

1610    SNUM  =  SNUM  +  1 

1620  GOTO  1570 

1630  CLOSE  #1 

164  0  LOCATE  20,5:  PRINT  " 
ii 

1650  LOCATE  21,8:  PRINT  "  ENTER  1  TO  USE  THESE  AOA 

SELECTIONS,  2  TO  ENTER  NEW  SELECTIONS" 

1660  LOCATE  22,33:  INPUT  "OR  3  TO  ABORT" ;AN4 

1670  COLOR  1 

1680  ON  AN4  GOTO  2150,  1740,  5180 

1690  ' 

1700  ' 

1710  '  AOA  SELECTION  -  IF  PREVIOUS  AOA  VALUES  ARE  NOT 

DESIRED 

1720  ' 

1730  ' 

1740  CLS:  LOCATE  3,26:  COLOR  4:  PRINT  "ANGLE  OF  INCIDENCE 

SELECTION":  COLOR  1 
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1750  LOCATE  8,22:  PRINT  "ENTER  THE  DESIRED  ANGLE  OF 

INCIDENCE" 

1760  LOCATE  10,25:  PRINT  "LIMITS  ARE  +18  TO  -18  DEGREES 

1770  LOCATE  12,19:  PRINT  "ENTER  AS  MANY  AOA'S  AS  YOU  WISH  TO 

MEASURE" 

1780  LOCATE  15,31:  PRINT  "ENTER  -99  TO  STOP" 

1790  LOCATE  20,22:  PRINT  "TOTAL  NUMBER  OF  SELECTIONS 

ENTERED : " 

1800  COLOR  4 

1810  SNUM  =  0 

1820  • 

183  0  ■  THE  VARIABLE  SNUM  STANDS  FOR  "SELECT  NUMBER" 

1840  ' 

1850  WHILE  SELECT  <>  -99 

1860      LOCATE  2  0,58:  PRINT  SNUM 

PRINT  "      " 

INPUT  "AOA:"; SELECT 

1 

=  SELECT 


1870 
1880 
1890 
1900 


LOCATE  16,43: 
LOCATE  16,37: 
SNUM  =  SNUM  + 
WANTAOA(SNUM) 


STORE  SELECTED  AOA  VALUES  FOR  FUTURE  USE 


1910  WEND 

1920  LOCATE  22,11:  INPUT  "ENTER  1  TO  STORE  THESE  FOR  LATER 

USE,  0  NOT  TO  STORE  THEM";AN6 

1930  IF  AN6  <>  0  GOTO  1990  ELSE  2150 

1940 

1950 

1960 

1970 

1980 

1990  OPEN  "AOA. DAT"  FOR  OUTPUT  AS  #1 

2000  SNUM  =  1 

2010  WHILE  WANTAOA(SNUM)  <>  -99 

2020    PRINT  #1,  USING  "###" ;WANTAOA(SNUM) 

2030    SNUM  =  SNUM  +  1 

2040  WEND 

2050  CLOSE  #1 

2060  LOCATE  22,11:  PRINT  " 


2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 


LOCATE  22,24:  PRINT  "ENTER  1  TO  CONTINUE,  2  TO  ABORT" 
LOCATE  23,28:  INPUT  "OR  3  TO  GO  BACK  ONE  PAGE";  AN7 
ON  AN7  GOTO  2150,  5180,  1470 


POSITION  TURNTABLE  TO  SELECTED  AOA 


EXPER  =  1 

WANTAOA  =  WANTAOA( EXPER) 

CLS:  COLOR  20:  LOCATE  10,35:  PRINT  "WAIT.  .  ." 

COLOR  1 

IF  EXPER  =  1  THEN  GOTO  22  00  ELSE  GOTO  224  0 
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2200  LOCATE  14,22:  PRINT  "AIRFOIL  IS  POSITIONING  TO  FIRST 

AOA" 

2210  GOSUB  6400 

2220  GOSUB  6530 

2230  GOTO  2270 

2240  LOCATE  14,23:  PRINT  "AIRFOIL  IS  POSITIONING  TO  NEXT 

AOA" 

2250  GOSUB  6400 

2260  GOSUB  6530 

2270  ' 

2280  ' 

2290  '  ADJUST  WIND  TUNNEL  VELOCITY  WITHIN  THE  CALIBRATED  . 

RANGE 

2300  ' 

2310  ' 

2320  GOTO  2340 

23  3  0  COLOR  1:  LOCATE  15,*  48:  PRINT  "UPDATING  VELOCITY" 

2340  GOSUB  6940 

2350  GOSUB  7310 

23  60  COLOR  1 

2370  CLS:  LOCATE  6,15:  PRINT  "ADJUST  TUNNEL  SPEED  WITHIN  THE 

CALIBRATED  RANGE  OF" 

2380  COLOR  4:  LOCATE  8,10:  PRINT  USING  "###.###" ;BNUM1 

2390  LOCATE  8,18:  PRINT  "FT/SEC":  COLOR  1:  LOCATE  8,25: 

PRINT  " ( " 

2400  LOCATE  8,26:  PRINT  USING  "###.###" ;BNUM3 

2410  LOCATE  8,33:  PRINT  "  KTS)  TO":  COLOR  4:  LOCATE  8,42 

242  0  PRINT  USING  "###.###" ;BNUM2 :  LOCATE  8,50:  PRINT 

"FT/SEC" 

2430  COLOR  1:  LOCATE  8,57:  PRINT  "(":  LOCATE  8,58:  PRINT 

USING  "###.###" ;BNUM4 

2440  LOCATE  8,66:  PRINT " KTS ) " 

2450  LOCATE  12,28-:  PRINT  "PRESENT  TUNNEL  VELOCITY:" 

2460  COLOR  4 

2470  LOCATE  15,41: 

2480  LOCATE  14,33: 

2490  LOCATE  14,41: 

2500  LOCATE  16,33: 

2510  LOCATE  16,41: 

2520  Ql  =  Q 

2530  LOCATE  20,23: 


PRINT  "  " 

PRINT  USING  "###.##.#"  ;VEL 

PRINT  "FT/SEC" 

PRINT  USING  "###.###", 'KTS 

PRINT  "KNOTS" 

PRINT  "ENTER  1  TO  UPDATE  TUNNEL  VELOCITY" 


2540  COLOR  1:  LOCATE  21,24:  PRINT  "(UPDATE  TAKES  APPROX  13 

SECONDS) " 

2550  COLOR  4:  LOCATE  23,32:  INPUT  "OR  2  TO  CONTINUE" ;AN10 

2560  IF  AN10  <>  2  GOTO  2330  ELSE  2590 

2570 

2580 

2590  '  ADJUST  POSITION  OF  THE  TRAVERSE  ASSEMBLY 

2600 

2610 

2  62  0  COLOR  1 
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2  63  0  NUMBER  =  0:  READING  =  0 
2640  WHILE  READING  <>  -99 

2650     IF  NUMBER  =  0  GOTO  2660  ELSE  2690 

2660     CLS:  LOCATE  6,16:  PRINT  "MOVE  THE  PITOT- 

STATIC/HOTWIRE  TRAVERSE  ASSEMBLY" 

2670     LOCATE  8,16:  PRINT  "TO  THE  POSITION  FOR 

YOUR  FIRST  READING" 

2680     COLOR  4:  LOCATE  8,23:  PRINT  "FOUR  INCH":  COLOR  1: 

GOTO  2710 

2690     CLS:  LOCATE  6,16:  PRINT  "MOVE  THE  PITOT- 

STATIC/HOTWIRE  TRAVERSE  ASSEMBLY" 

2700     LOCATE  8,29:  PRINT  "FOR  YOUR  NEXT  READING" 

2710     LOCATE  12,27:  PRINT  "PRESENT  POSITION:        INCH" 

2720     INPUT. CH  =  4 

2730     CALL  OUTPUT (RELAY. MUX. 01, INPUT. CH) 

2740     CALL  ENABLE. OUTPUT (RELAY. MUX. 01) 

2750     CALL  MEASURE (DMM. 01, INCH) 

2760     POSITION (NUMBER)  =  (INCH*1000)/1. 053 

2770     COLOR  4 

2780     LOCATE  12,45:  PRINT  USING  "##.##"; POSITION (NUMBER) 

2790     LOCATE  20,27:  PRINT  "PRESS  ANY  KEY  WHEN  YOU  ARE" 

2800     LOCATE  21,22:  PRINT  "THROUGH  ADJUSTING  TRAVERSE 

POSITION" 

2810     COLOR  1 

2820     A$  =  INKEY$:  IF  A$  =• ""  THEN  GOTO  2750  ELSE  2830 

2830     CLS:  COLOR  20:  LOCATE  10,36:  PRINT  "WAIT...":  COLOR 

1 

2840     LOCATE  14,19:  PRINT  "HOT  WIRE  VOLTAGE  DATA  NOW  BEING 

COLLECTED" 

2850 

2860 

2870 

2880 

2890 

2900     INPUT. CH  =  1 

2910     CALL  OUTPUT (RELAY. MUX. 01, INPUT. CH) 

2920     1=0:  SUM1  =  0 

2930     FOR  I  =  1  TO  40 

2940  .      IF  I  =  1  GOTO  2950  ELSE  2970 

2950        TIME  =  .5 

2960        CALL  DELAY (TIME) 

2970        CALL  MEASURE (DMM. 01, G (I) ) 

2980        SUM1  =  SUM1  +  ABS(G(I)) 

2990     NEXT  I 

3  000     AVREAD1  =  SUM1/4  0 
3  010     HWVOLT (NUMBER)  =  AVREAD1 
3020     CORRFAC  =  (TS  -  TCAL)/(TS  -  TF) 
3030     ' 
3040     ' 

3  050     •  COLLECT  PITOT-STATIC  SYSTEM  DATA 
3060     • 
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COLLECT  HOTWIRE  VOLTAGE  DATA 


3070     ' 

3080     GOSUB  6940 

3  090     GOSUB  7310 

3100     '  PS.VEL  IS  THE  TUNNEL  VELOCITY  AS  MEASURED  BY  THE 

PITOT-STATIC  TUBE 

3110     '  S.PR  IS  THE  STATIC  PRESSURE  TAKEN  FROM  THE  STATIC 

PRESSURE  PORT 


3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

1.03 

3250 

3260 

1.015 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

1.015 

3370 


PS.VEL (NUMBER)  =  VEL 

i 


GET  VELOCITY  PROFILE  FROM  HOTWIRE  DATA 


X (NUMBER)  =  POSITION (NUMBER) 
ECORR  =  (HWVOLT (NUMBER) A2)  *  CORRFAC 
V2 (NUMBER)  =  ((ECORR  -  B)  /  A)A2 
VI  =  V2(0) 

Y (NUMBER)  =  V2 (NUMBER) /Vl 
IF  WANTAOA  >=  10  GOTO  3240  ELSE  3260 
IF  Y (NUMBER)  <  .97  GOTO  3270  ELSE  IF  Y (NUMBER)  > 
GOTO  3270  ELSE  Y (NUMBER)  =  1 
GOTO  3270 

IF  Y (NUMBER)  <  .985  GOTO 
GOTO  3270  ELSE  Y (NUMBER)  = 
YINT (NUMBER)  =  Y (NUMBER) 
IF  ABS(YINT( NUMBER) )  <= 


3270 

1 

*  (1  - 

01  THEN 


ELSE  IF  Y (NUMBER)  > 

Y (NUMBER))   " 
YINT (NUMBER)  =0 


GET  VELOCITY  PROFILE  FROM  PITOT-STATIC  SYSTEM  DATA 


VINF  =  PS.VEL(O) 

Y2 (NUMBER)  =  PS .VEL (NUMBER) /VINF 

IF  Y2 (NUMBER)  <  .985  GOTO  3370  ELSE  IF  Y2 (NUMBER) 
GOTO  3370  ELSE  Y2 (NUMBER)  =  1 

Y2INT (NUMBER)  =  (PS .VEL (NUMBER ) /VINF) - 
( ( PS. VEL (NUMBER)/ VINF) A2) 
3380     IF  ABS(Y2INT (NUMBER) )  <= 
3390     ' 
3400     ' 

3410     CLS:  BEEP:  COLOR  4:  LOCATE  11,28: 
TO  CONTINUE" 

3420     LOCATE  13,14:  INPUT  "OR  ENTER  -99 
COMPLETED  YOUR  LAST  READING" /READING 


01  THEN  Y2 INT (NUMBER)  =  0 


PRINT  "PRESS  ENTER 


IF  YOU  HAVE 


3430 
3440 
3450 
3460 
3470 
3480 


IF  READING  =  -99  GOTO  3450  ELSE  3440 
NUMBER  =  NUMBER  +  1 
COLOR  1 


WEND 

EXPER  =  EXPER  +  1 

CLS:  LOCATE  10,15 


PRINT  "DATA  VALUES  ARE  BEING  STORED 


IN  DATA  FILES  LABELED" 
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3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 

3660 

3670 

3680 

3690 

3700 

V2/VINF 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

DISPL 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

V2/VINF 

3920 

3930 

3940 

3950 

3960 


COLOR  4:  LOCATE  12,35:  PRINT  "PSTAT.DAT" 


LOCATE  14,35:  PRINT  "HWIRE.DAT" 

COLOR  1 

OPEN  "PSTAT.DAT"  FOR  OUTPUT  AS  #1 

J  =  0 

FOR  J  =  1  TO  NUMBER 

PRINT  #1,USING"###.###";X(J),  Y2 ( J) 
NEXT  J 
CLOSE  #1 

OPEN  "HWIRE.DAT"  FOR  OUTPUT  AS  #1 
J  =  0 
FOR  J  =  1  TO  NUMBER 

PRINT  #1,USING"###.###";X(J) ,  Y(J) 
NEXT  J 
CLOSE  #1 
COLOR  1:  CLS 
LOCATE  2,22 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 


2,53 
3,42 
4.31 


LOCATE  1,31:  PRINT  "EXPERIMENTAL  DATA" 
PRINT '"ANGLE  OF  INCIDENCE  (DEGREES)  =" 
PRINT  USING  "###.#"; WANT AOA 
PRINT  "VELOCITY  PROFILES" 
PRINT  "HOTWIRE" 
PRINT  "PITOT-STATIC" 
PRINT  "POSITION  (IN) 
V2/VINF" 


4,56 
5,7: 

V2  (FT/SEC) 
J  =  0:  COLOR  4 
FOR  J  =  1  TO  NUMBER 

IF  J  >=  30  GOTO  4080  ELSE  3740 
IF  J  =  29  GOTO  3990  ELSE  3750 
IF  J  >=  16  GOTO  3920  ELSE  3760 
IF  J  =  15  GOTO  3830  ELSE  3770 


V2  (FT/SEC) 


LOCATE  J+6,11: 
LOCATE  J+6,25: 
LOCATE  J+6,40: 
LOCATE  J+6,53: 
LOCATE  J+6,68: 
GOTO  4140 
COLOR  1:  LOCATE 
MORE  DATA" 


PRINT  USING  »##.##»;X(J) 
PRINT  USING  "###.###»;V2(J) 
PRINT  USING  »#.####» ;Y(J) 
PRINT  USING  "###.###"; PS. VEL( J) 
PRINT  USING  "#.###";Y2(J) 


22,23:  PRINT  "PRESS  ANY  KEY  TO 


A$  =  INKEY$: 
COLOR  1:  CLS 
LOCATE  2,22 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
V2 


2,53 
3,42 
4,31 
4,56 
5,7: 


IF  A$  =  ""  THEN  3840 

LOCATE  1,31:  PRINT  "EXPERIMENTAL  DATA" 
PRINT  "ANGLE  OF  INCIDENCE  (DEGREES)  =" 
PRINT  USING  "###.#" ;WANTAOA 
PRINT  "VELOCITY  PROFILES" 
PRINT  "HOTWIRE" 
PRINT  "PITOT-STATIC" 
PRINT  "POSITION  (IN)    V2  (FT/SEC) 


(FT/SEC) 


COLOR  4 

LOCATE  J-8,11:  PRINT 

LOCATE  J-8,25:  PRINT 

LOCATE  J-8,40:  PRINT 

LOCATE  J-8,53:  PRINT 


V2/VINF" 

USING  »##.##";X(J) 

USING  »###.###»;V2(J) 

USING  »#.####»;Y(J) 

USING  "###.###";PS.VEL(J) 


86 


3970    LOCATE  J-8 , 68 :  PRINT  USING  "#.####" ;Y2 (J) 
3980    GOTO  4140 

3990    COLOR  1:  LOCATE  22,23:  PRINT  "PRESS  ANY  KEY  TO 
DISPLAY  MORE  DATA" 


4000  A$  =  INKEY$: 

4010  COLOR  1:  CLS 

4  020  LOCATE  2 

403  0  LOCATE 

4040  LOCATE 

4050  LOCATE 

4060  LOCATE 

407  0  LOCATE 


22: 

53: 
42: 
31: 
56: 


V2/VINF 


5,7: 


IF  A$  =  ""  THEN  4000 
:  LOCATE  1,31:  PRINT  "EXPERIMENTAL  DATA" 

PRINT  "ANGLE  OF  INCIDENCE  (DEGREES)  =" 

PRINT  USING  "###.#" ;WANTAOA 

PRINT  "VELOCITY  PROFILES" 

PRINT  "HOTWIRE" 

PRINT  "PITOT-STATIC" 
PRINT  "POSITION  (IN)    V2  (FT/SEC) 


4080 
4090 
4100 
4110 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
COPY 
4210 


V2  (FT/SEC) 


COLOR  4 

LOCATE  J-22,11 
LOCATE  J-22,25 
LOCATE  J-2  2,4  0 
LOCATE  J-22,53 
LOCATE  J- 2 2, 68 
NEXT  J 


V2/VINF" 


PRINT  USING  "##.##";X(J) 
PRINT  USING  »###.###";V2(J) 
PRINT  USING  "#.####»;Y(J) 
PRINT  USING  "###.###" ;PS.VEL(J) 
PRINT  USING  "#.####";Y2(J) 


SEND  DATA  SUMMARY  TO  LASER  PRINTER 


PRINT  "WOULD  YOU  LIKE  A  HARD 


INPUT  "ENTER  1  FOR  YES,  2  FOR 


COLOR  1:  LOCATE  22,20: 

OF  THIS  DATA?" 

COLOR  4:  LOCATE  23,27: 
N0";AN11 
4220  IF  AN11  <>  2  GOTO  4230  ELSE  4570 

CLS:  LOCATE  12,26:  PRINT  "DATA  SENT  TO  LASER  PRINTER" 

COLOR  1:  TIME$  =  "00:00:00" 

TIME  =  VAL(RIGHT$(TIME$,2) ) 

IF  TIME  <  5  GOTO  4250  ELSE  4270 

LPRINT 

LPRINT  TAB (31):  LPRINT  "EXPERIMENTAL  DATA" 

LPRINT 

LPRINT  TAB (22):  LPRINT  "ANGLE  OF  INCIDENCE  (DEGREES) 


4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 


4310  LPRINT  TAB(53); 

4320  LPRINT  USING  "###.#" ;HAVEAOA 

43  30  LPRINT 

434  0  LPRINT  TAB (42) :  LPRINT 

4350  LPRINT  TAB (31):  LPRINT 

STATIC" 

43  60  LPRINT  TAB (7):  LPRINT 

V2/VINF      V2  (FT/SEC) 

4370  LPRINT 

4380  J  =  0 

4390  FOR  J  =  0  TO  NUMBER 

4400    .LPRINT  TAB (11) ; 


"VELOCITY 
"HOTWIRE 

"POSITION 
V2/VINF" 


PROFILES" 


PITOT- 


(IN)   *  V2  (FT/SEC) 
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4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 


LPRINT  USING"##.##" ;X(J) ; 

LPRINT  TAB (25) ; 

LPRINT  USING"###.###";V2(J) ; 

LPRINT  TAB (40) ; 

LPRINT  USING"#.###";Y(J) ; 

LPRINT  TAB (53) ; 

LPRINT  USING" ###.###"; PS. VEL (J) ; 

LPRINT  TAB (68) ; 

LPRINT  USING"#.###";Y2(J) 
NEXT  J 
LPRINT  CHR$(12) 


PRINT  VIEW  PLOT  OF  HOTWIRE  VELOCITY  PROFILE 


COLOR  1 
CLS:  LOCATE 


12,9:  PRINT  "WOULD  YOU  LIKE  TO  SEE  A  PLOT 


OF  THE  HOTWIRE  VELOCITY  PROFILE?" 

4590  COLOR  4:  LOCATE  14,27:  INPUT  "ENTER  1  FOR  YES,  2  FOR 

NO";AN12 

4600  COLOR  1:  IF  AN12  <>  2  GOTO  4610  ELSE  4640 

CLS:  SHELL" PART21.BAT" 
i 


4610 

4620 

4630 

4640 

PRINTER 

4650  • 

4660  ' 

4670  CLS:  LOCATE 

THIS  GRAPH?" 

4680  LOCATE 

NO";AN13 

4690  COLOR  1:  IF 

4700  GOTO  4830 

4710 

4720 

4730 

4740 

4750 

4760  CLS:  COLOR  1:  LOCATE  10,21:  PRINT  "THIS  PROGRAM  WILL 

AUTOMATICALLY  LEAVE" 

4770  LOCATE  12,22:  PRINT  "BASICA  AND  ENTER  A  GRAPHICS 

PROGRAM . " 

4780  LOCATE  14,14:  PRINT  "THE  PLOT  TAKES  APPROXIMATELY  4-6 


SEND  HARD  COPY  OF  HOTWIRE  VELOCITY  PROFILE  TO  LASER 


12,19:  PRINT  "WOULD  YOU  LIKE  A  HARD  COPY  OF 
14,27:  COLOR  4:  INPUT  "ENTER  1  FOR  YES,  2  FOR 
AN13  <>  2  THEN  GOSUB  4730  ELSE  GOTO  4860 


SUBROUTINE  FOR  INSTRUCTIONS  ON  LASER  PRINTER 


COMPLETE" 
4:  LOCATE 


21,27:  PRINT  "PRESS  ANY  KEY  TO 


MINUTES  TO 

4790  COLOR 

CONTINUE" 

4800  A$  =  INKEY$:  IF  A$  =  ""  THEN  4800 

4810  COLOR  1 

482  0  RETURN 
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PRINT  VIEW  PLOT  OF  PITOT-STATIC  VELOCITY  PROFILE 


483  0  CLS:  SHELL  "PART22.BAT" 

4840  ■ 

4850  ' 

4860  ' 

4870  ' 

4880  ' 

4890  CLS:  LOCATE  12,7:  PRINT  "WOULD  YOU  LIKE  TO  SEE  A  PLOT 

OF  THE  PITOT-STATIC  VELOCITY  PROFILE?" 

4900  COLOR  4:  LOCATE  14,27:  INPUT  "ENTER  1  FOR  YES,  2  FOR 

NO";AN14 

4910  COLOR  1:  IF  AN14  <>  2  GOTO  4920  ELSE  4950 

4920  CLS:  SHELL  "PART2  3.BAT" 

4930  ' 

4940  ' 

4950  '  SEND  HARD  COPY  OF  PITOT-STATIC  VELOCITY  PROFILE  TO 

LASER  PRINTER 

4960  ' 

4970  ' 

4980  CLS:  LOCATE  12,19:  PRINT  "WOULD  YOU  LIKE  A  HARD  COPY  OF 

THIS  GRAPH?" 

4990  LOCATE  14,27:  COLOR  4:  INPUT  "ENTER  1  FOR  YES,  2  FOR 

NO";AN15 

5000  COLOR  1:  IF  AN15  <>  2  THEN  GOSUB  4730  ELSE  GOTO  5040 

5010  CLS:  SHELL  "PART24.BAT" 

5020  ! 

5030  ' 

5040  '  PRINT  CONCLUDING  SCREEN 

5050  ■ 

5060  ' 

5070  CLS:  COLOR  4:  LOCATE  6,12:  PRINT 
I'********************************************************" 

5080  LOCATE  10,12:  PRINT 


5090  LOCATE  7,12:  PRINT  "*" 
5100  LOCATE  9,12:  PRINT  "*" 
5110  LOCATE  8,67:  PRINT  "*" 


LOCATE  8,12:  PRINT  "*" 
LOCATE  7,67:  PRINT  "*" 
LOCATE  9,67:  PRINT  "*" 

5120  COLOR  1:  LOCATE  8,16:  PRINT  "THIS  CONCLUDES  THE  HOT 

WIRE  FLOW  MEASUREMENT  LAB" 

5130  COLOR  4:  LOCATE  20,18:  PRINT  "ENTER  1  TO  REPEAT  THIS 

LAB  WITH  THE  NEXT  AOA" 

5140  LOCATE  21,15:  PRINT  "ENTER  2  TO  CONTINUE  WITH  DRAG 

COEFFICIENT  ANALYSIS" 

5150  LOCATE  22,20:  INPUT  "ENTER  3  TO  EXIT  BASICA  AND  RETURN 

TO  DOS";AN18 

5160  ON  AN18  GOTO  2160,  5190,  5180 

5170  GOTO  5190 

5180   CLS:  SYSTEM 

5190  XI  =  X(l):  X2  =  X(NUMBER-3) 

5200  1=0 

5210  FOR  I  =  0  TO  NUMBER 

5220    Y(I)  =  YINT(I) 
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523  0  NEXT  I 
5240  GOSUB  8380 
5250  PROFILE1  =  Z 
5260  CHORD  =  12 
5270  • 

5280  '  CDU1  IS  THE  UNCORRECTED  DRAG  COEFFICIENT  FOR  THE 
HOTWIRE  SYSTEM 
5290  ' 

5300  CDU1  =  ((2/CHORD)*PROFILEl) 
5310  »' 

5320  '  SBF  IS  THE  SOLID  BLOCKING  FACTOR 
5330  ' 

5340  SBF  =  .00525 
5350  ' 

53  60  '  WBF  IS  THE  WAKE  BLOCKAGE  FACTOR 
5370  WBF  =  .13333  *  CDU1 
5380  ' 

5390  •  CD1  IS  THE  CORRECTED  DRAG  COEFFICIENT  FOR  THE  HOTWIRE 
DATA  ANALYSIS 
5400  ' 

5410  CD1  =  CDU1*(1-(3*SBF)-(2*WBF) ) 
5420  ' 

543  0  '  FIND  DRAG  COEFFICIENT  FOR  THE  PITOT-STATIC  SYSTEM 
5440  ' 
5500  1=0 

5510  FOR  I  =  0  TO  NUMBER 
5520    Y(I)  =  0 
553  0  NEXT  I 
5540  1=0 

5550  FOR  I  =  0  TO  NUMBER 
5560    Y(I)  =  Y2INT(I) 
5570  NEXT  I 

5580  IF  Y (NUMBER)  <  .01  THEN  Y (NUMBER)  =  0 
5590  GOSUB  8380 
5600  PROFILES  =  Z 
5610  ■ 

5620  '  CDU2  IS  THE  UNCORRECTED  DRAG  COEFFICIENT  FOR  THE 
PITOT-STATIC  SYSTEM 
5630  ' 

5640  CDU2  =  ((2/CHORD)*PROFILE3) 
5650  WBF  =  .13333*CDU2 
5660  ' 

5670  '  CD2  IS  THE  CORRECTED  DRAG  COEFFICIENT  FOR  THE  PITOT- 
STATIC  SYSTEM 
5680  • 

5690  CD2  =  CDU2*(1-(3*SBF)-(2*WBF) ) 
5700  ' 
5710  ■ 

5720  '  PRINT  DRAG  COEFFICIENT  SUMMARY  ON  SCREEN 
5730  ' 
5740  ' 
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5750  COLOR  1 
SUMMARY" 
5760  LOCATE 
5770  LOCATE 
5780  COLOR  1 
5790  LOCATE 
5800  LOCATE 
5810  COLOR  4 
5820  LOCATE 
5830  COLOR  1 
5840  LOCATE 
5850  LOCATE 
5860  COLOR  4 
5870  LOCATE 
5880  LOCATE 
Cd  SUMMARY?" 
5890  LOCATE 
5900  IF  AN19 
5910 
5920 
5930 
5940 
5950 
5960  LPRINT 
5970  LPRINT 
5980  LPRINT 
5990  LPRINT 

6000  LPRINT 

—  ii  • 

—  i 

6010  LPRINT 
6020  LPRINT 
603  0  LPRINT 
6040  LPRINT 
6050  LPRINT 
6060  LPRINT 
6070  LPRINT 
6080  LPRINT 
6090  LPRINT 
6100  LPRINT 
6110  LPRINT 
6120  LPRINT 

613  0  LPRINT 

614  0  LPRINT 
6150  LPRINT 
6160  LPRINT 
6170  LPRlNT 
6180  LPRINT 
6190  LPRINT 
62  00  LPRINT 
6210  LPRINT 
622  0  LPRINT 


:  CLS:  LOCATE  2,28:  PRINT  "DRAG  COEFFICIENT 

4,26:  PRINT  "AIRFOIL  AT        DEGREES  AOA" 

4,37:  COLOR  4:  PRINT  USING  "###.#" ;HAVEAOA 

:  LOCATE  7,33:  PRINT  "HOTWIRE  SYSTEM" 

9,28:  PRINT  "UNCORRECTED  Cd  =" 

11,29:  PRINT  "CORRECTED  Cd  =" 

:  LOCATE  9,45:  PRINT  USING  "#.####" ;CDU1 

11,44:  PRINT  USING  "#.####" ;CD1 

:  LOCATE  14,30:  PRINT  "PITOT-STATIC  SYSTEM" 

16,28:  PRINT  "UNCORRECTED  Cd  =" 

18,29:  PRINT  "CORRECTED  Cd  =" 

:  LOCATE  16,45:  PRINT  USING  "#.####" ;CDU2 

18,44:  PRINT  USING  "#.####";CD2 

21,17:  PRINT  "WOULD  YOU  LIKE  A  HARD  COPY  OF  THE 

22,27:  INPUT  "ENTER  1  FOR  YES,  2  FOR  NO" ;AN19 
<>  2  GOTO  5930  ELSE  6290 


SEND  DRAG  SUMMARY  TO  LASER  PRINTER 


TAB(28):  LPRINT  "DRAG  COEFFICIENT  SUMMARY" 


TAB (22) :  LPRINT  "ANGLE  OF  INCIDENCE  (DEGREES) 

TAB (53) ; 

US ING  "###.#»; HAVE AOA 


TAB (33):  LPRINT  "HOTWIRE  SYSTEM" 

TAB(28) :  LPRINT  "UNCORRECTED  Cd  =" ; 

TAB (45) ; 

USING  "#.####" ;CDU1 

TAB (29):  LPRINT  "CORRECTED  Cd  =" ; 

TAB (44) ; 

USING  "#.####";CD1 


TAB(30):  LPRINT  "PITOT-STATIC  SYSTEM" 

TAB (28) : LPRINT  "UNCORRECTED  Cd  =" ; 

TAB (45) ; 

USING  "#.####"; CDU2 
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6230 
6240 
6250 
6260 
6270 
6280 
6290 


LPRINT  TAB(29):  LPRINT  "CORRECTED  Cd  =" 

LPRINT  TAB (44) ; 

LPRINT  USING  "#.####" ;CD2 

LPRINT  CHR$(12) 


LOCATE  21,16:  PRINT  " 


63  00  LOCATE  22,25:  PRINT  " 
ii 

6310  LOCATE  21,18: 

THE  NEXT  AOA" 

6320  LOCATE  22,20: 

TO  DOS";AN20 

6330  IF  AN20  <>  2 

6340  STOP 

6350 

6360 

6370 

6380 

6390 

6400 

6410 

6420 

6430 

6440 


PRINT  "ENTER  1  TO  REPEAT  THIS  LAB  WITH 
INPUT  "ENTER  2  TO  EXIT  BASICA  AND  RETURN 


GOTO  1430  ELSE  5180 


SUBROUTINE  FOR  READING  AOA 


INPUT. CH  =  5 

CALL  OUTPUT (RE LAY. MUX. 01, INPUT. CH) 

CALL  ENABLE . OUTPUT ( RELAY. MUX. 01) 

TIME  =  .5 

CALL  DELAY (TIME) 
6450  CALL  MEASURE (DMM. 01, AOA) 
6460  IF  AOA  <  0  GOTO  6470  ELSE  6490 

HAVEAOA  =  (AOA-.003)*1000/10 

GOTO  6500 

HAVEAOA  =  (AOA-.003)*1000/10 

RETURN 


6470 
6480 
6490 
6500 
6510 
6520 
6530 
6540 
6550 
6560 
6570 
6580 
6590 
6600 
6610 
6620 
6630 
6640 


10)  /  10)  >=  .25 


SUBROUTINE  FOR  AOA  SELECTION 


WHILE  (CINT(ABS(WANTAOA  -  HAVEAOA)  ' 

ABSDIFF  =  ABS(WANTAOA  -  HAVEAOA) 

IF  ABSDIFF  >  36  THEN  ABSDIFF  =36 

SWITCH. 1  =  1 

SWITCH. 2  =  2 

CALL  OPEN . CHANNEL ( RELAY . ACT .01, SWITCH . 1 ) 

CALL  OPEN. CHANNEL (RE LAY. ACT. 01, SWITCH. 2) 

IF  WANTAOA  >  HAVEAOA  GOTO  6670  ELSE  6800 
i 

•  THIS  PART  OF  THE  SUBROUTINE  MOVES  THE  TURNTABLE 


6650 

FORWARD 

6660    ' 

6670    IF  (CINT (ABSDIFF  * 

ELSE  TIMING  =  1 


10)  /  10)  >=  .4  THEN  TIMING  =  100 
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6680  COUNT  =  ( (ABSDIFF  *  TIMING)  +  190)  /  1000 

6690  IF  COUNT  >  3  THEN  COUNT  =  3 

6700  CALL  CLOSE. CHANNEL (RE LAY. ACT. 01, SWITCH. 1) 

6710  CALL  DELAY (COUNT) 

6720  CALL  OPEN. CHANNEL (RELAY. ACT. 01, SWITCH. 1) 

6730  TIME  =  .5 

6740  CALL  DELAY (TIME) 

6750  GOSUB  6400 

6760  GOTO  6900 

6770  ' 

6780  '  THIS  PART  OF  THE  SUBROUTINE  MOVES  THE  TURNTABLE  IN 

REVERSE 

6790  ' 

6800  IF  (CINT (ABSDIFF  *  10)  /  10)  >=  .4  THEN  TIMING  =  100 

ELSE  TIMING  =  1 

6810  COUNT  =  ((ABSDIFF  *  TIMING)  +  190)  /  1000 

6820  IF  COUNT  >  3  THEN  COUNT  =  3 

683  0  CALL  CLOSE. CHANNEL (RELAY. ACT. 01, SWITCH. 2) 

6840  CALL  DELAY (COUNT) 

6850  CALL  OPEN. CHANNEL (RELAY. ACT. 01, SWITCH. 2) 

6860  TIME  =  .5 

6870  CALL  DELAY (TIME) 

6880  GOSUB  6400 

6890  GOTO  6900 

6900  WEND 

6910  RETURN 

6920 

6930 

6940 

6950 

6960 

6970  COLOR  1 

6980  Al  =  38.709457# 

6990  A2  =  .037085566# 

7000  A3  =  5.649552E-05 

7010  BO  =  .10086091# 

7020  Bl  =  25727. 94369# 

7030  B2  =  -767345. 8295# 

7040  B3  =  78025595. 81# 

7050  B4  =  -9247486589# 

7060  B5  =  6.97688E+11 

7070  B6  =  -2.66192E+13 

7080  B7  =  3.94078E+14 

7090  B8  =  0 

7100  B9  =  0 

7110  REF.OUT  =  8 

7120  CALL  OUTPUT (RELAY. MUX. 01, REF.OUT) 

7130  TIME  =  .5 

714  0  CALL  DELAY (TIME) 

7150  CALL  MEASURE (DMM. 01, V) 

7160  TR  =  100  *  V 
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SUBROUTINE  FOR  MEASURING  WIND  TUNNEL  TEMPERATURE 


ER  =  (TR*(A1+TR*(A2+TR*A3) ) )*10^-6 

CALL  SET.RANGE(DMM.01,R2  00MILLI) 

TC  =  7 

CALL  OUTPUT (RELAY. MUX. 01 ,TC) 

TIME  =  .5 

CALL  DELAY (TIME) 

CALL  MEASURE (DMM. 01, ET) 

E  =  ER  +  ET 

Z  =  B5+E*(B6+E*(B7+E*(B8+E*B9) ) ) 

T  =  B0+E*(B1+E*(B2+E*(B3+E*(B4+E*Z) ) ) ) 

TF  =  (1.8  *  T)  +32 

RETURN 


7170 

7180 

7190 

7200 

7210 

7220 

7230 

7240 

7250 

7260 

7270 

7280 

7290 

7300 

7310 

7320 

7330 

7340 

7350 

7360 

7370 

7380  CLS:  COLOR  1: 

THE  PARASCIENTIFIC 

7390 

7400 


SUBROUTINE 


FOR  WIND  TUNNEL  DATA: 
STATIC  PRESS  IN  TEST  SECTION 
AIR  DENSITY  IN  TEST  SECTION 
DYNAMIC  PRESSURE  IN  TEST  SECTION 
TEST  SECTION  VELOCITY,  FT/SEC  AND 


KTS 


LOCATE  10,15: 
TRANSDUCER" 
COLOR  4:  LOCATE  10,30:  PRINT  "P2" 
LOCATE  12,15:  PRINT  "MAKE  SURE  NO 


PRINT  "DEPRESS  BUTTON 


ON 


COLOR  1 
OTHER  FUNCTION 


PRINT 
""  THEN 


BUTTONS  ARE  DEPRESSED" 

7410  COLOR  4:  LOCATE  21,28: 

7420  A$  =  INKEY$:  IF  A$  = 

7430  1=0:  SUM3  =  0 

7440  FOR  I  =  1  TO  10 

7450    CALL  MEASURE (DIG. IN. 01, J (I) ) 

7460    SUM3  =  SUM3  +  J(I) 

7470  NEXT  I 

7480  DN  =  SUM3/10 

7490  GOSUB  7900 

7500  PSTAT  =  PRESS  *  144 

7510  IF  PSTAT  <  1000  THEN  RETURN 

7520  TTEMP  =  TF  +  459.67 

7530  GASCONST  =53.3 

7540  DENSITY  =  PSTAT/ (GASCONST  * 

7550  CLS:  COLOR  1:  LOCATE  10,15: 

THE  PARASCIENTIFIC  TRANSDUCER" 


"PRESS 
7420 


ENTER  TO  CONTINUE" 


TTEMP) 

PRINT  "DEPRESS 


BUTTON 


ON 


COLOR  1 
OTHER  FUNCTION 


7560  COLOR  4:  LOCATE  10,30:  PRINT  "PI": 

7570  LOCATE  12,15:  PRINT  "MAKE  SURE  NO 

BUTTONS  ARE  DEPRESSED" 

7580  COLOR  4:  LOCATE  21,28:  PRINT  "PRESS  ENTER  TO  CONTINUE" 

7590  A$  =  INKEY$:IF  A$  =  ""  THEN  7590 

7600  1=0:  SUM1  =  0 

7610  FOR  I  =  1  TO  70 

7620    CALL  MEASURE (DIG. IN. 01, F (I) ) 

7630    F(I)  =  F(I)  +  .007 
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7640    SUM1  =  SUM1  +  F(I) 

7650  NEXT  I 

7660  DN  =  SUM1/70 

7670  GOSUB  7900 

7680  VELOC  =  PRESS 

7690  IF  VELOC  >  10  THEN  VELOC  =  0 

7700  DELTAP  =  VELOC* 14 4 

7710  Q  =  DELTAP/. 93 

7720  IF  Q  <  1  THEN  Q  =  0 

7730  GC  =  32.174 

7740  VEL  =  SQR(ABS((2*Q*GC)/DENSITY) ) 

7750  KTS  =  VEL  *  .5924 

7760  FILE$  =  "HOTWIRE. HPC" 

7770  CALL  INITIALIZE. SYSTEM (FILE$) 

7780  CALL  ENABLE . SYSTEM 

7790  CALL  SET. FUNCTION (DMM. 01, DCVOLTS) 

78  00  CALL  SET.RANGE(DMM:01,AUTOM) 

7810  CALL  SET.SPEED(DMM.01,R2.5) 

7820  CALL  DISABLE. INT. TRIGGER (DMM. 01) 

783  0  CALL  ENABLE. OUTPUT (RELAY. MUX. 01) 

7840  CALL  SET. NUM.  BITS  (DIG.  IN.  01,R16) 

7850  LEVEL  =  5 

7860  CALL  SET. THRESHOLD (DIG. IN. 01, LEVEL) 

7870  RETURN 

7880  ■ 

7890  ' 

7900  ■  SUBROUTINE  FOR  CONVERTING  DECIMAL  PRESSURE  VALUES  TO 

NUMERIC  VALUES 

7910  ' 

7920  ■ 

7930  COUNT  =  0:  PRESS  =  0 

7940  ' 

7950  '  DN  IS  THE  MEASURED  DECIMAL  NUMBER 

7960  '  DEC  IS  THE  DECIMAL  TRANSFORM  MATRIX 

7970  '  PRESS  IS  THE  CONVERTED  NUMERIC  PRESSURE  VALUE 

7980  ' 

7990  DEC(16)  =  10:  DEC(15)  =  8:  DEC(14)  =  4:  DEC(13)  =  2: 

DEC (12)  =  1 

8000  DEC(ll)  =  .8:  DEC(IO)  =  .4:  DEC(9")  =  .2:  DEC(8)  =  .1: 

DEC(7)  =  .08 

8010  DEC(6)  =  .04:  DEC(5)  =  .02:  DEC(4)  =.  .01:  DEC(3)  = 

8.000001E-03 

8020  DEC(2)  =  .004:  DEC(l)  =  .002 

8030  ' 

8  04  0  IF  DN  <  0  THEN  PRESS  =  DEC (16)  ELSE  GOTO  8060 

8050  COUNT  =  -32768! 

8060  IF  (DN-COUNT)/16384  >=  1  THEN  PRESS  =  PRESS  +  DEC (15) 

ELSE  GOTO  8080 

8070  COUNT  =  COUNT  +  16384 

8080  IF  (DN-COUNT)/8192  >=  1  THEN  PRESS  =  PRESS  +  DEC (14) 

ELSE  GOTO  8100 
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8090  COUNT  =  COUNT  +  8192 

8100  IF  (DN-COUNT)/4096  >= 

ELSE  GOTO  812  0 

8110  COUNT  =  COUNT  +4096 

8120  IF  (DN-COUNT)/2048  >= 

ELSE  GOTO  8140 

8130  COUNT  =  COUNT  +  2048 

8140  IF  (DN-COUNT)/1024  >= 

ELSE  GOTO  8160 

8150  COUNT  =  COUNT  +  1024 

8160  IF  (DN-COUNTJ/512  >= 

ELSE  GOTO  8180 

8170  COUNT  =  COUNT  +512 

8180  IF  (DN-COUNT)/256  >= 

GOTO  8200 

8190  COUNT  =  COUNT  +  256 

8200  IF  (DN-COUNT)/128  :>= 

GOTO  8220 

8210  COUNT  =  COUNT  +128 

8220  IF  (DN-COUNT)/64  >= 

GOTO  8240 

823  0  COUNT  =  COUNT  +  64 

8240  IF  (DN-COUNT)/32  >= 

GOTO  8260 

8250  COUNT  =  COUNT  +  32 

8260  IF  (DN-COUNT)/16  >= 

GOTO  8280 

8270  COUNT  =  COUNT  +  16 

8280  IF  (DN-COUNT)/8  >= 

GOTO  8300 

8290  COUNT  =  COUNT  +  8 

8300  IF  (DN-COUNT)/4  >= 

GOTO  8320 

8310  COUNT  =  COUNT  +  4 

8320  IF  (DN-COUNT)/2  >= 

GOTO  8340 

83  30  COUNT  =  COUNT  +  2 

8340  IF  (DN-COUNT)  = 

8350  RETURN 

8360 

8370 

8380 

8390 


1  THEN  PRESS  =  PRESS  +  DEC (13) 


1  THEN  PRESS  =  PRESS  +  DEC (12) 


1  THEN  PRESS  =  PRESS  +  DEC (11) 


1  THEN  PRESS  =  PRESS  +  DEC (10) 


1  THEN  PRESS  =  PRESS  +  DEC (9)  ELSE 


1  THEN  PRESS  =  PRESS  +  DEC (8)  ELSE 


1  THEN  PRESS  =  PRESS  +  DEC (7)  ELSE 


1  THEN  PRESS  =  PRESS  +  DEC (6)  ELSE 


1  THEN  PRESS  =  PRESS  +  DEC (5)  ELSE 


1  THEN  PRESS  =  PRESS  +  DEC (4)  ELSE 


1  THEN  PRESS  =  PRESS  +  DEC (3)  ELSE 


1  THEN  PRESS  =  PRESS  +  DEC (2)  ELSE 


1  THEN  PRESS  =  PRESS  +  DEC(l) 


SUBROUTINE  FOR  INTEGRATION 
INPUTS:    X(I), 


Y(I)   TABLE  VALUES  IN 


11  PART  2.  DAT" 

8400 

8410 

8420 

8430 

8440 

8450  V  =  NUMBER 

8460  GOSUB  9060 


V  =  NUMBER  OF  TABLE  VALUES 
XI,  X2   INTEGRATION  RANGE 
Z  =  CALCULATED  INTEGRAL 
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8490  RETURN 

8500 

************************************************************ 

***** 

8510  ■  AKIMA  SPLINE  FITTING  SUBROUTINE 

8520  ■    THE  INPUT  TABLE  IS  (X(I),  Y(I)),  WHERE  Y(I)  IS  THE 

DEPENDENT  VARIABLE 

8530  '    THE  INTERPOLATION  POINT  IS  X,  WHICH  IS  ASSUMED  TO 

BE  IN  THE  RANGE  OF 

8540  '    THE  TABLE  WITH  AT  LEAST  ONE  TABLE  POINT  TO  THE 

LEFT,  AND  THREE  TO  THE 

8550  '    RIGHT.  Y  IS  RETURNED  AS  THE  INTERPOLATED  VALUE.  N 

IS  RETURNED  AS  AN 

8560  '    ERROR  CHECK  (N  =  0  IMPLIES  AN  ERROR) .  DIMENSION 

M,X,Y,  AND  Z  IN  THE 

8570  '    CALLING  PROGRAM. 

8580  N  =  1 

8590  '  CHECK  TO  SEE  IF  X  IS  IN  THE  TABLE  RANGE 

8600  IF  X  >  =  X(l)  GOTO  8630 

8610  N  =  0 

8620  RETURN 

8630  IF  X  <=  X(V-3)  GOTO  8660 

8640  N  =  0 

8  650  RETURN 

8660  X(0)  =  2*X(1)-X(2) 

8  670  '  CALCULATE  AKIMA  COEFFICIENTS 

8680  1=0 

8690  FOR  I  =  1  TO  V-l 

8700   '  SHIFT  I  TO  1+2 

8710    M(I+2)  =  (Y(I+1)-Y(I) )/(X(I+l)-X(I)  ) 

8720  NEXT  I 

8730  M(V+2)  =  2*M(V+1)  -  M(V) 

8740  M(V+3)  =  2*M(V+2)  -  M(V+1) 

8750  M(2)  =  2*M(3)  -  M(4) 

8760  M(l)  =  2*M(2)  -  M(3) 

8770  1=0 

8780  FOR  I  =  1  TO  V 

8790    A  =  ABS(M(I+3)  -  M(I+2)) 

8800    B  =  ABS(M(I+1)  -  M(I) ) 

8810    IF  A+B  <>  0  GOTO  8840 

8820    Z(I)  =  (M(I+2)  +M(I+l))/2 

8830    GOTO  8850 

8840    Z(I)  =  (A*M(I+1)  +  B*M(I+2))/(A+B) 

8850  NEXT  I 

8860  ■  FIND  RELEVANT  TABLE  INTERVAL 

8870    1=0 

8880    1=1+1 

8890  IF  X  >=  X(I)  GOTO  8880 

8900  1=1-1 

8910  '  BEGIN  INTERPOLATION 

8920  B  =  X(I+1)  -  X(I) 
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8930  A  =  X  -  X(I) 

8940  Y  =  Y(I)+Z(I)*A+(3*M(I+2)-2*Z(I)-Z(I+l) ) *A*A/B 

8950  Y  =  Y+(Z(I)+Z(I+l)-2*M(I+2) )*A*A*A/(B*B) 

8960  RETURN 

8970  ' 

********************************************************** 

**** 

GENERAL  INTEGRATION  SUBROUTINE  (ITEG) 
INTERPOLATION  BY  AKIMA 

INTEGRATION  BY  ENHANCED  TRAPEZOIDAL  RULE 
THE  INTEGRATION  RANGE  IS  (XI,  X2) . 


8980 

8990 

9000 

9010 

9020  • 

LEAST 

9030  • 

X2. 

9040  ' 

9050  ' 

ERROR. 

9060 

9070 

9080 

9090 

9100 

9110 

9120 

9130 

9140 

9150 

9160 

9170 

9180 

9190 

9200 

9210 

9220 

9230 

9240 

9250 

9260 

9270 

9280 

9290 

9300 

9310 

9320 

9330 

9340 

9350 


ONE 


IT  IS 
TABLE 
VALUE 


ASSUMED  THAT  XI  <  X2 ,  AND  THAT  THERE  IS  AT 
TO  THE  LEFT  OF  XI,  AND  THREE  TO  THE  RIGHT  OF 


THE  RESULT  IS  RETURNED  IN  Z. 
AN  ERROR  CHECK  IS  RETURNED  IN 


Zl.  Zl  =  0  IMPLIES 


Z  =  0 
Zl  =  0 
1  CHECK 
IF  XI  < 
IF  X2  > 
1  IF  XI 
IF  XI  < 
=  XI 
=  X2 
=  X3 
=  1 
X2  = 
START 


TO  SEE  IF  END  POINTS  ARE  IN  ALLOWABLE  RANGE 

X(l)  THEN  RETURN 

X(V-3)  THEN  RETURN 

>  X2  THEN  SWITCH  AND  SET  FLAG 

X2  GOTO  9170 


X3 

XI 

X2 

Zl 

IF 
i 


FIRST 


XI  THEN  RETURN 
TRAPEZOIDAL  INTEGRATIONS 
INTEGRATION  TO  GET  II 

GOSUB  9380 

1  SECOND  ROUND  TO  GET  12 

GOSUB  9680 

1  RICHARDSON  EXTRAPOLATION 

Z  =  4*12/3  -  11/3 

•  CHECK  TO  SEE  IF  THE  END  POINTS 

IF  Zl  =  0  GOTO  9310 

Z  =  -Z 

X2  =  XI 

XI  =  X3 

1  RESET 

Zl  =  1 

RETURN 

1  Z  IS 

RETURN 
i 


HAVE  BEEN  REVERSED 


ERROR  FLAG 


THE  INTEGRAL  DESIRED 


************************************************************ 
* 

93  60  '  ROUTINE  FOR  THE  FIRST  TRAPEZOIDAL  INTEGRATION,  II 
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9370  '   Nl  KEEPS  TRACK  OF  THE  NUMBER  OF  INTERVALS 

9380  II  =  0 

9390  Nl  =  0 

9400  X  =  XI 

9410  '  FIND  THE  BEGINING  OF  THE  INTERVAL 

9420  '  GO  TO  BRANCH  WHICH  CALLS  THE  INTERPOLATION  ROUTINE 

9430  •  FIND  THE  INTERVAL,  I,  AND  THE  LEFT  END  POINT,  Y. 

9440  GOSUB  9990 

9450  IF  X2  >  X(I+1)  GOTO  9540 

9460  Nl  =  Nl  +  1 

9470  D  =  Y 

9480  X  =  X2 

9490  '  FIND  END  POINT  Y  VALUE 

9500  GOSUB  9990 

9510  II  =  (Y+D)*(X2-Xl)/2 

952  0  RETURN 

9530  '  AT  LEAST  ONE  TABLE  INTERVAL  MUST  BE  SUMMED  OVER 

9540  Jl  =  I 

9550  II  =  II  +  (Y+Y(I+1) )*(X(I+l)-X)/2 

9560  IF  X2  <  X(Jl+3)  GOTO  9610 

9570  Nl  =  Nl  +  1 

9580  II  =  Il+(Y(Jl+l)+Y(Jl+3) )*(X(Jl+3)-X(Jl+l) )/2 

9590  Jl  =  Jl  +  2 

9600  GOTO  9560 

9610  X  =  X2 

9620  GOSUB  9990 

9630  II  =  I1+(Y+Y(J1+1) )*(X2-X(J1+1) )/2 

9640  Nl  =  Nl  +  1 

9650  RETURN 

9660  ' 

*********************************************************** 

*** 

9670  '  INTEGRATION  FOR  12 

9680  12  =  0 

9690  X  =  XI 

9700  GOSUB  9990 

9710  D  =  Y 

9720  IF  X2  >  X(I+1)  GOTO  9810 

9730  X  =  Xl+(X2-Xl)/2 

9740  GOSUB  9990 

9750  12  =  I2+(D+Y)*(X2-Xl)/4 

9760  D  =  Y 

9770  X  =  X2 

9780  GOSUB  9990 

9790  12  =  I2+(D+Y)*(X2-Xl)/4 

98  00  RETURN 

9810  X  =  Xl+(X(I+l)-Xl)/2 

9820  Jl  =  I 

9830  GOSUB  9990 

9840  12  =  I2+(Y+D)*(X-Xl)/2 

9850  12  =  I2+(Y+Y(J1+1) ) * (X ( Jl+1) -X) /2 
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9860  IF  X2  <  X(Jl+2)  GOTO  9900 

9870  12  =  I2+(Y(Jl+l)+Y(Jl+2) ) * (X( Jl+2) -X( Jl+1) )/2 

9880  Jl  =  Jl  +  1 

9890  GOTO  9860 

9900  X  =  X2-(X2-X(J1+1) )/2 

9910  GOSUB  9990 

9920  D  =  Y 

9930  12  =  I2+(Y(Jl+l)+D)*(X2-X)/2 

9940  X  =  X2 

9950  GOSUB  9990 

9960  12  =  I2+(D+Y)*(X2-X(J1+1) )/4 

9970  RETURN 

9980  '  GO  TO  AKIMA  SPLINE  INTERPOLATION  SUBROUTINE 

9990  GOSUB  8580 

10000  '  JUST  RETURNED  FROM  INTERPOLATION  SUBROUTINE 

10010  '  RETURN  TO  PROGRAM 

10020  RETURN 

10030  END 
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